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I. INTRODUCTION 


LL modern theories of reaction velocity 

recognize the importance of precise de- 
terminations of the energy of activation. For this 
purpose the temperature coefficients of the 
velocity of gaseous reactions and of non-ionic 
reactions in solution have been studied in some 
detail. 

Comparable data for ionic reactions are ex- 
ceedingly meager and the theoretical implications 
have been virtually neglected. The activation 
energies of ionic reactions are of particular 
interest for the following reasons. 

The partial molal heat capacities of ions are 
large negative! 2 values, which depend upon the 
electric charge, the total concentration, and 
certain specific factors arising from interaction 
with the solvent. Since the temperature de- 
pendence of the energy of activation involves 
heat capacities, it is conceivable for reactions 
between ions that Fact. may not be independent 
of temperature. The approximation that Eact. is 
constant is almost universally accepted in inte- 
grating the Arrhenius equation. 

Interionic attraction modifies the rates of ionic 
reactions by changing the number and character 
of the effective collisions, a circumstance which 


, ‘La Mer and Cowperthwaite, J. Am. Chem. Soc. 55, 1004 
1933), 


* Randall and Rossini, ibid. 51, 323 (1929). 


emphasizes the ‘necessity for considering the 
orientation of the molecules at collision.*: 4 

The experimental evidence that activation 
energies of ionic reactions are dependent upon 
both temperature and concentration will be given 
in later papers of this series. The present paper 
will treat the simpler case of non-ionic reactions 
from the viewpoint of the statistical foundations 
of the Arrhenius equation; the derivation and 
limitations of the Brénsted equation; the validity 
of collision theories; and the experimental evi- 
dence for the temperature dependence of Ea¢¢. in 
non-ionic reactions. 


II. THE ARRHENIUS EQUATION AND 
Its DERIVATIVES 


We start with the Arrhenius equation in the 


Bact. 
(1) 


which R. C. Tolman® has shown follows from 
statistical mechanics for a bimolecular reaction. 
To emphasize the salient features for later 
discussion, we give the simpler development valid 
for a unimolecular reaction and refer the reader 


3La Mer and Kamner, J. Am. Chem. Soc. 53, 2832 
(1931). 

4La Mer, Chem. Rev. 10, 179 (1932). 

5R. C. Tolman, Statistical Mechanics, Chem. Catalog 
Co., pp. 259-69, 1927. 
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to Tolman and Kassel’s monographs’ for details 
of the origin of the term 37 in the bimolecular 
case. 

_Consider a system of molecules of index 7, 
distributed in various quantum levels having the 
energies ¢; and the statistical weights p;. The 
molecules in the higher reactive quantum states 
will be indicated by the index j, and will be 
assigned undetermined specific individual rate 
(reactivity) constants k;. The number of mole- 
cules, ;, in the j*® quantum state is 


(2) 


where kT is the Boltzmann constant times the 
absolute temperature. 
The overall, measurable, specific rate constant 
ky, is then 
1 dn _ dik ipje 
k= —- (3) 
n dt > 


Logarithmic differentiation yields 


The differentiation yields: 
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dT Dik Dipie~ ei! 
= (1/kT?)[€—€ ] = Enct./ RT”. (5) 


The last term in the bracket is the average 
energy € of all molecules, while the first term refers 
to the average energy of the molecules in j states, in 
which the usual weight factor p; has been replaced 
by k;p;. This means that € refers to the average 
energy of the molecules which react. The quantity 
N(€—€) will be defined as the energy of activation 
per gram mole. 

Since the reacting molecules constitute but a 
minute fraction of the total, no sensible error is 
made by identifying € as the average energy of 
those molecules which do not react. 

To obtain an expression for the variation of 
the energy of activation with temperature, we 
differentiate Tolman’s expression, Eq. (4), and 
obtain a quantity Cact. which we shall call the 
heat capacity of activation. 


d[kT2(d In (6) 


Dd spies! ( De || 


dE = (1/kT?) (8) 


In the case of heat capacities, the bar refers to 
the partial molal quantity, whereas with the 
energies ¢ it refers to the average. 

The quantity [e?—é*] is the difference between 
the mean of the squares of the energies of all the 
molecules and the square of their mean energy, 
and consequently represents the fluctuations’: § 
in their energy. [e2—€*]/kT? is equal to the 
calorimetrically measurable partial heat capacity 
of all the molecules, which heat capacity is 
indistinguishable from that of the non-reacting 
molecules. 
By the same token [e?—e? ]/kT? represents the 


6L. S. Kassel, Kinetics of Homogeneous Gas Reactions, 
Chem. Catalog Co., pp. 19-26, 1932. 

7R. H. Fowler, Statistical Mechanics, Cambridge, 1929; 
Eq. (1435), p. 513. 

8 J. W. Gibbs, Statistical Mechanics, Longmans, 1928; 
(Chap. VII, Eq. (205), Chap. XIV and XV). 


(Cjract.) —C)) = Crct. (8a) 


partial heat capacity of the molecules which react, 
but differs from the thermodynamic heat capac- 
ity of the molecules in activated states since each 
p; has been replaced by kjp; in the averaging 
process. 

The integration of Eq. (1) under the assump- 
tion that E,c+. is independent of T is equivalent 
to setting Eq. (8) equal to zero, i.e., SCact. 
= Cx act.) —(Ca+Czg) is zero for the bimolecular 
reaction A+BeX OR. 

When heat capacities are neglected free energy 
and total energy become indistinguishable. This 
approximate treatment, corresponding in the 
integrated form to 


In k=2.3B—Eact./RT (9) 


has met with considerable success in the treat- 
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reat- 


ment of gaseous reactions, notably in the hands 
of Hinshelwood.° 

The correction which we introduce by con- 
sidering the temperature dependence of F,.¢. can 
be understood better by poiriting out the close 
analogy to the corresponding thermodynamic 
process of calculation. The integral 


Lact. Bact. OE act. 
f dT=— + f ( a In 7. 
RT oT 


The last term by definition is equal to Cact.dln T, 
which is in complete analogy with the 
thermodynamic expression for the entropy S; 
namely, S= {(dE/dT)d in T. We therefore con- 
sider it entirely justifiable to speak of the 
entropy of activation, S.ct., provided it is kept 
clearly in mind that S,.:. contains the kinetic 
factors k; found in the first terms of Eq. (7). 

It cannot be emphasized too strongly, that to 
neglect the entropy of activation, as is the case, 
(Eq. 9), when E£,-¢. is taken independent of 7, 
has no more justification than would be the 
omission of the entropy terms involved in the 
calculation of a thermodynamic equilibrium. 

In the case of electrolytes the assumption that 
Cuct.=0 is certainly not justified a priori. We 
shall accordingly integrate under the assumption 
that E,.:. is a function of 7, even though we 
cannot take full advantage of the improvement 
produced; thus, 


Fact. 
In k=constant+3 In (10) 
RT 


The free energy of activation, Fact. =Eact. 
—TSxet.= Frxeact. —(Fat Fr), where the bars 
refer to the chemical potentials. The subscript, 
activation, is retained for the X or critical state 
to emphasize that Fx act.) contains the k; factors 
and therefore refers to the molecules which 
react, whereas F, and Fy, are free from k; 


°C. N. Hinshelwood, Kinetics of Homogeneous Gas 
Reactions. Oxford, 1929, p. 42, where the view is advanced 
that “if when log & is plotted against 1/T a straight line is 
not obtained, this is an almost certain indication that the 
observed reaction is a composite one made up of two or 
More concurrent reactions differently influenced by 
temperature.” Cf. Kassel, reference 6, pp. 149-56, where 
the data on hydrogen iodide have been subjected to a 
critical reexamination and the conclusion drawn that 
Exct, increases with temperature. 
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factors. The X state is the rate determining in- 
termediate state. 

The constant in (10) is evaluated by com- 
parison with In k when F,,;.=0. For a process 
without activation, every collision will be 
fruitful and the constant=In Z°—} In T. 
For a dilute gas, Z° is the gas kinetic constant 
for the number of collisions per unit of time 
= RT(m\+mz2)/m mz }', and has the di- 
mensions of a frequency (time). 


III. SysTEMS WITH INTERACTION. THE 
BRONSTED EQUATION 


To account for changes in the medium with 
changes in concentration, Bronsted generalized 
the classical theory of reaction velocity by 
introducing the activity coefficients of the re- 
actants and. the critical complex in terms of a 
kinetic factor fafe/fx. His theory although most 
useful in treating ionic reactions is by no means 
so limited. The foundation of his equation has 
been the subject of much discussion! and 
a direct derivation is highly desirable. 

The statistics, upon which Eqs. (1-7) are 
founded, were developed primarily for perfect 
gas mixtures. For such systems, the quantities 
€;, €; and k; are independent of both temperature 
and concentration. However, when we consider 
more concentrated systems, and particularly the 
liquid state, ¢«;, ¢; and presumably also &; will 
depend upon the mutual energy of the molecules. 
These quantities may consequently depend upon 
the concentration and temperature, since the 
energy of each molecular system depends upon 
the configuration of all other molecular systems 
in the assembly. 

The ensemble" method of Gibbs, in contrast to 
the classical development of Boltzmann recog- 
nizes the fundamental character of the problem 
of interaction energy by including it in the 
conditional equations employed in deriving the 


10 Bronsted, Zeits. f. physik. Chemie 102, 169 (1922); 
115, 337 (1925). 

1 Bjerrum, Zeits. f. physik. Chemie 108, 82 (1924); 118, 
251 (1925). 

12 Christiansen, Zeits. f. physik. Chemie 113, 35 (1924). 

18 Skrabal, Zeits. f. physik. Chemie B3, 247 (1929). 

4 W. Gibbs, reference 8, pp. 36-7. Consult also Fowler, 
reference 7, where certain aspects of the problem are 
stated but not solved, pp. 164-5, 181-2, 307. 
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distribution laws. Unfortunately an _ explicit 
solution by the Gibbs method encounters un- 
surmountable mathematical difficulties, so that 
for the present we must abandon a rigorous 
statistical derivation and resort to thermo- 
dynamic analogies as our guide. In a thermo- 
dynamic procedure the correction factors due to 
interaction will enter as the empirically de- 
termined activity coefficients. 

In this section we shall drop, temporarily, the 
subscript, activation, to simplify the notation, 
but without losing sight of the kinetic character 
of such quantities. In analogy with thermody- 
namic procedure we decompose F into an ideal 
part Fig, and a non-ideal or excess part, F, 
arising from interaction effects. 

Thus 


F=F(T)+RT In f(T, c) (11) 
=F. + (11a) 


equations which serve as definitions for the 
activity coefficient, f. The small c will refer to 
concentration, the large capital C to ideal heat 
capacity. 

The omission of the term RT In c for Fia. 
requires justification. In evaluating the free 
energy of activation, we are concerned only with 
the individual energy levels which are inde- 
pendent of concentration for the perfect gas 
mixture. Terms of the type RT In c arise from 
the disorder which occurs when molecules of 
different types are mixed but not when like 
types are mixed (Gibbs paradox). In the present 
instance, the integration producing Fia, involves 
only quantities which are characteristic of the 
energies of each individual group of molecules, 
and consequently no terms of the type RT In ¢ 
can appear in the expression for F. This may be 
more evident by considering the formal analogy 
to the previously mentioned thermodynamic 
integral S= {Cd |In T, where C represents the 
ideal heat capacity, which of course cannot 
produce the term R In c, which arises from the 
Gibbs paradox. 

On the other hand, the RT In f term originates 
from the integration of the corresponding non- 
ideal heat capacity term produced by interaction 
effects whereby an excess entropy and a corre- 
sponding excess free energy result. 
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Introducing (11) and (11a) in (10) yields 
Ink=In (12) 
—[Fxcact.) —(Fat Fa) ]\ 


13 
(13) 


k=Z° exp 


As a result, we have Bronsted’s Eq. (14) 

k=Ro( fafe)/fx (14) 
which can also be written 

k= koe 
when & is set equal to 
—[Fx°— ] 

Z° exp ( 

RT 


The choice of the constants 
which are independent of concentration, implies 
that fafs/fx is normalized as unity for the 
reference solvent; hence any dependence upon 
concentration arising from a change in the 
medium (for example, neutral salt addition) 
finds expression in terms of the activity coeff- 
cients. In passing from the gaseous to new liquid 
states of reference, any purely kinetic factors such 
as viscosity, which can affect the numerical value 
of Z° will have to be taken into account.’ 
ky thus corresponds to the specific velocity 
constant in the reference solvent when freed 
from all interactions between A, B and X except 
idealized collisions, but not from interaction with 
the solvent. If relative activity coefficients are 
normalized for each new solvent, the effect of 
interaction with the solvent will appear in the 
values of Fx°, F4°, Fg° required for each solvent. 

We believe that in our method of approach we 
have given for the first time a clear meaning to 
the quantities fy or Fx°, by empaAasizing the 
limitations which must be placed upon all 
attempts to interpret fx as a pure thermodynamic 
quantity, for example, when fx is replaced by a 
simple theoretical expression like the Debye 
limiting law which obviously cannot contain f; 
factors. 

Only in the special case where the k; values are 
all equal or only one j level is involved, do the &; 
values cancel in the averaging process, so that fr 


4 Jowett, Phil. Mag. 8, 1059 (1929). ’ 
16 Moelwyn-Hughes, J. Chem. Soc. (London) 1932, 9: 


) (15) 
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can be interpreted as a pure thermodynamic 
quantity. Equality of the k; factors is unlikely. 
On the other hand, if the reaction proceeds 
through a single quantum state j, then by Eqs. 
(7) and (8), the heat capacity of the reacting 
molecules is zero and Cact.=0; i.e., Eact. must be 
constant or decrease with rising temperature. 
Consequently if Eact. increases with T we have 
convincing evidence that the reaction proceeds 
through more than one reactive level and the X 
state is multiple. 

As a corollary of our theory we have the 
proposition that whenever £,,;¢. rises with tem- 
perature, substitution of a thermodynamic ex- 
pression (f. ex. (D.H.L.L.)) for fx is illegitimate. 
Even if Cact.=0, we cannot be certain that the X 
state is not multiple, since a finite heat capacity 
for the reacting molecules may be counter- 
balanced by the heat capacity of the non- 
reacting molecules. 

Hinshelwood’s statement® that inconstancy of 
Exct. indicates a composite reaction is, of course, 
true when considered from the standpoint of 
paths involving different quantum levels, but 
from the context his statement is obviously 
intended for the appearance of side reactions 
involving the production of new chemical species 
or entirely new mechanisms in the chemical 
sense. 


IV. Entropy oF ACTIVATION, AND THE 
COLLISION THEORY 


Recently Moelwyn-Hughes" has extended the 
simple collision theory of gaseous reactions to 
solutions. He determines Exact. from measure- 
ments of the temperature coefficient, and from 
the gas kinetic value of Z, using reasonable values 
of the molecular radii o:, and o2 he computes a 
theoretical value of k. 

For a large number of reactions the calculated 
value of k differs from the observed by a factor of 
only 10 to 100. Consequently Moelwyn-Hughes 
concludes that the simple collision theory is 
applicable to solutions. 

There are two reasons why this is not con- 
Vincing: (1) For a group of reactions (mostly 

* Moelwyn-Hughes, Chem: Rev. 10, 241 (1932); Nature 
129, 316 (1932); J. Chem. Soc. (1932) 95; with Hinshel- 
Wood, (1932) 230; with Rolfe, (1932) 241. 
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bimolecular associations), Moelwyn-Hughes finds 
k is 104 to 10° times smaller than that calculated, 
while our unpublished results for the 6-bromo- 
propionate-thiosulfate reaction yield rates which 
are many times greater than those calculated 
from his theory, in spite of the fact that the 
repulsive action of the electric charges might be 
expected to reduce the rate. (2) Moelwyn- 
Hughes’ theory neglects completely the possi- 
bility of an entropy of activation (Sact.) resulting 
from a temperature dependence of Eact.. 
By comparing Eqs. (9) and (10), we find 


B=log 29+ Sact./2.3R, (16) 


where Sact. is equal to 


Cuct.d In T. 
‘ 0 


The Hinshelwood-Moelwyn-Hughes theory 
makes the tacit assumption that the heat 
capacities of activation of the reacting and of the 
non-reacting molecules remain identical down to 
absolute zero. There is no valid theoretical reason 


‘for such an assumption, especially when one 


considers that for such minor differences as that 
of proton spin, the heat capacities of ortho and 
para hydrogen differ considerably at low temper- 
atures,'® and consequently possess different 
entropies at room temperature. 

There are six important cases connecting 
| 2° and B. 

Case I. Cact. >0 at all temperatures; E,.¢. rises 
with increasing 7; Sact. >0, and B>log Z°. 

Case II. Cact.>0 in the range of experimental 
measurement (high 7) but Cacr.KO at lower 
temperatures. Then Sact.<0 and B<log Z° in 
the range of investigation. 

Case III. Cact. <0 at all temperatures. Eact. 
decreases; Sact. <0; B<In Z°. 

Case IV. Cact.<0 at high 7, but Cact.>>0 at 
low T. Sat.>0 and B>log Z°® in range of 
measurement. 

Case V. Condition for the validity of Arrhenius 
assumption. Cet. =0 within experimental error in 


temperature range of measurement; Sact.>0 or 


18 Eucken, Hiller and Clusius, Zeits. f. physik. Chemie 
B4, 142 (1929); Treatise on Physical Chemistry, Van 
Nostrand (1931), p. 305. 
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S<0 depending upon behavior of Cact. at low T. TABLE I. The decomposition of acetone-dicarboxylicacid 


B>log Z°, or B<log Z°. in water.* 

Case VI. Condition for the validity of the (1st order); time units minutes™. 
Hinshelwood-Moelwyn-Hughes theory. Cact.=0 joe B Enct Cot 
at all temperatures or varies in such a way at low oa ae poe ; 
temperatures that =0; then B=log 2°. 4 

10-20 $4.77 22993 | 
35 
V. EXPERIMENTAL DATA 20-30 14.03 23340 | 
49 

It is unfortunate that E,ct. cannot be meas- 30-40 14.37 23826 = 

ured at low temperatures to reveal towhatextent 409-50 14.13 23477 | ‘ 
—23 
the neglect of the entropy of activation accounts P se ones 


for the failure of the collision theory. However, a 
critical examination of the best existing data 
discloses that Enact. varies with T in the range of 
measurement. All values of Exact. have been 
recalculated from original data by using R 
= 1.9885 calories and T = 273.13 for the ice point. 
B isin minutes. 

Table I refers to Wiig’s results for the de- 
composition of acetone-dicarboxylic acid from 0 
to 60°C at 10° intervals. No statement of the 
reproducibility of the individual measurements is 
given, but since with two different analytical 
techniques k checks to 1 percent (+300 calories 
in Exct.) the probable error in a series for the 
same technique may be less. Dr. Wiig informs me 
that a reexamination of his original data shows 
that the average deviation from the mean at each 
temperature is about 1 percent. 

The reversal in sign of Cact. would be disre- 
garded were it not uniformly exhibited by 
Segaller’s extensive studies on the reaction be- 
tween primary, iso- and secondary alkyl iodides 
with sodium phenoxide. Our calculations for the 


Solvent, absolute ethanol. B(2nd order); B=min.—. 


TABLE II. The reaction of sodium phenoxide with primary alkyl iodides.* 


* Wiig, J. Phys. Chem. 32, 961 (1928); 34, 596 (1930). 


primary iodides, which cover a 50° range are 
given in Table II. Segaller called attention to the 
fact that his measurements at higher temperature 
did not fall on a straight line for a log k, (1/7) 
plot, and contends the discrepancy cannot be due 
to experimental error. 

In Fig. 1 are plotted the values of Cac. 
calculated from Wiig’s data (J. Am. Chem. Soc. 
53, 4729 (1930) Table VI) for kes°/kis° and 
k35°/ko5° for the decomposition of citric acid in 
concentrated sulfuric acid. With the exception of 
the two points at low water content, it is very 
interesting to note that with increasing water 
content, Cact. changes quite uniformly from 
negative to positive values. This means that the 
addition of water increases the heat capacity of 
the reactive complex. 

In Table III are given calculations from G. M. 
Murphy’s measurements on the diacetone- 
alcohol reaction, considered as a bimolecular 


Alkyl t=30.1°-42.5° t=42.5°-58.5° t=58.5°-80.1° 
iodide E B E B E B 
Methyl 20,945 13.28 21,876 13.92 (931) — — ae 
i Ethyl 21,067 12.71 . 21,604 13.08 537 21,263 12.86 —341 
Propyl 21/508 12.60 22/048 —«:12.98 540 21561 12.65 —487 
| Butyl 20,560 11.92 22,197 13.05 637 21,551 12.62 — 646 
Amyl 20,696 11.63 22,353 12.79 657 21,385 12.15 —968 
Hexyl 20,622 «11.94 22'175 ~—-13.01 553 21,770 12.75 —405 
z Heptyl 20,834 12.07 22,183 13.00 549 21,958 12.85 —225 
Octy 21,199 12.31 22'364 13.12 (1,165) 22/050 —314 
4 Hexadecyl 21,307 12.37 22,188 12.98 (881) 22,146 12.96 —42 
Average 41 cals. —2to ~= 
Cact. cals. deg. 


* D. Segaller, J. Chem. Soc. 105, 106 (1914). 
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Fic. 1. Heat capacity of activation in calories/deg. mole 
(15-25°; 25-35°) for the water inhibited reaction. Thermal 
decomposition of citric acid in concentrated sulfuric acid 
as function of water content. Data by Wiig, J. Am. Chem. 
Soc. 53, 4729 (1930), Table VI. 


TaBLE III. The decomposition of diacetone-alcohol by 0.1M 
sodium hydroxide. E precise to + 0.7 percent.* 


B(mols. liter~! min.-' is computed on the basis of a 
bimolecular reaction between NaOH and diacetone- 
alcohol. B(ist order)=B(2nd order)—1. Log Z2°(¢ 
=5A°) = 13.36 (min.~). 


Solvent E B E B 

Water+0.1M NaOH = 20°-30° 25°-35° Cact. 
0% CH;OH 17,040 12.19 18,098 12.94 211 
10% “ 17,638 12.50 18,670 13.20 206 
20% 18,311 12.85 19,481 13.70 234 
30% 18,525 12.88 19,694 13.74 234 
0% “ 18,852 12.99 19,822 13.71 194 
50% 19,349 13.23 20,080 13.78 146 
600% 19,776 13.42 20,578 14.03 160 


*G. M. Murphy, J. Am. Chem. Soc. 53, 977 (1931). 


reaction with 0.1M NaOH. A precision of 0.7 
percent or 126 calories is claimed for Enct.. It 
will be noted that E,.:. increases over a ten 
degree interval by roughly 1000 calories or eight 
times the probable error. A calculation based 
upon 5° intervals indicates decreasing values of 
Cot, with increasing J as in Table II. This 
reaction is now being investigated in this 
laboratory over wider ranges of temperature to 
verify this deduction. 


VI. ConcLUDING REMARKS 


After this paper was practically completed for 
Publication the interesting papers of Scheffer" 


a and Brandsma, Rec. Trav. Chim. 45, 522 
6). 
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and Brandsma”’ were discovered, where the im- 
portance of considering the free energy and 
entropy of the intermediate state is emphasized. 
Their viewpoint differs from ours in at least one 
important respect. They consider that these pure 
thermodynamic properties determine the rate, 
whereas we consider the rate determining prop- 
erties of the intermediate state to be in general 
kinetic quantities, since they contain k; factors. 
Except for this difference, their thesis that part of 
the entropy of activation arises from ‘“‘steric 
factors” and ‘“‘sensitive districts” in the molecule 
is one which we had developed independently 
from our studies of the temperature dependence 
of ionic reactions, considered from the standpoint 
of the orientation theory.”! 

We can also support Scheffer and Brandsma’s 
conclusion that the effect of a catalyst does not 
always result in a lowering of the energy of 
activation, by the following example. The ve- 
locity of the sodium bromoacetate-sodium thio- 
sulfate reaction in highly dilute solution is 
increased over sixfold by the addition of small 
amounts of lanthanum chloride.” 

For the range 0—25°C, it has been found that 
Fact. increases from 15,850+64 to 17,700+55 
calories by the addition of 0.03M LaCl;. An 
increase of 1850 calories in E,-¢. should of itself 
decrease the velocity, were it not more than 
counterbalanced by an increase in the entropy 
of activation as exhibited by an increase in B 
from 11.20 to 13.40. As was anticipated quite 
different results, are found for the 6-bromo- 
propionate-thiosulfate reaction. 

The author wishes to acknowledge the con- 
structive criticisms which he has received in the 
preparation of this paper from his colleagues and 
from Professor Otto Halpern of the Department 
of Physics of New York University. 


VII. ConcLusions 


(1) Evidence is presented that the energy of 
activation, measured by the temperature coeffi- 
cient of the reaction velocity, varies with temper- 


20 Brandsma, ibid. 47, 94 (1928); 48, 1205 (1929). 

"1 Reference 3, p. 2850; reference 4, p. 210. 

22La Mer and Fessenden, J. Am. Chem. Soc. 54, 2351 
(1932). 

23 Unpublished work of La Mer and Kamner. 
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ature. This necessitates a consideration of the 
activation quantities, heat capacity, free energy 
and entropy. 

(2) The kinetic character of these quantities 
becomes evident when they are developed from 
R. C. Tolman’s statistical formulation of the 
Arrhenius equation. 

(3) The Brénsted equation follows readily 
from the concept of free energy of activation. 
~However, the development shows that the ac- 
tivity coefficient of the critical complex fx;act.) is 
in general a kinetic quantity which becomes 
equivalent to the corresponding thermodynamic 
quantity only when all quantum levels have the 
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same reactivity constant, or when only one 
quantum level is involved. 

(4) A positive heat capacity of activation 
furnishes evidence that the critical state is 
multiple in the quantum sense and that it is 
illegitimate to replace fx(act.. by a thermo- 
dynamic expression. 

(5) The Hinshelwood-Moelwyn-Hughes col- 
lision theory can be valid only when the entropy 
of activation is zero, which is improbable from 
the standpoint of quantum theory, the thermo- 
dynamic requirements of equilibrium processes, 
or the available data for activation energies in 
the ranges of experimental measurement. 
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Evaluation of the Series Which Arise in the Calculation of Thermodynamic Quantities 
From Spectroscopic Data 
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(Received February 13, 1933) 


In cases where knowledge of the molecular spectrum of a 
gas is sufficiently detailed to include the “stretching” and 
“interaction” terms in the expression for the rotational- 
vibrational energy states, calculation of the free energy, 
entropy and heat capacity for high temperatures becomes 
very laborious when the usual approximations are em- 
ployed. In this paper, methods are developed which reduce 


many such calculations to a matter of interpolation in 
appended tables; for cases which lie outside the range of 
the tables, approximations are given which materially 
reduce the labor of computation. The procedure is illus- 
trated by calculations for chlorine, bromine, hydrogen 
chloride, carbon monoxide, oxygen and nitric oxide. 


FE calculation of thermodynamic quanti- 
ties (heat capacities, entropies and ‘‘free 
energies’) from spectroscopic data is now a 
familiar process and provides an easy and 
accurate method of predicting chemical equi- 
libria. In carrying out such calculations for 
diatomic molecules, one of two alternative pro- 
cedures has in general been adopted. The first of 
these consists in computing the “state sum”’ and 
its derivatives from the empirically determined 
energy levels as presented in Eq. (3) below, the 
contribution of each rotational-vibrational state 
being calculated and the separate contributions 
added to form the sum. This method was first 
used by Hicks and Mitchell,! and has been 
extensively employed by Giauque and his associ- 
ates. It possesses the immense advantage of 
avoiding all assumptions as to the distribution of 
the levels, but becomes impossibly laborious if 
certain of the spectroscopic constants are small or 
the temperature is high. The second procedure is 
to treat the rotational-vibrational levels as those 
of a rigid rotator combined with those of either a 
harmonic or anharmonic oscillator, i.e., the terms 
involving D and F in Eq. (3) are ignored and the 
constant B is assumed to be independent of the 
vibrational quantum number. With these simpli- 
fying assumptions the rotational and vibrational 
parts of the state sum can be treated separately 


‘Hicks and Mitchell, J. Am. Chem. Soc. 48, 1520 (1926). 


Phim example—Giauque, J. Am. Chem. Soc. 52, 4816 
930). 


and the summation over the rotational levels can 
be replaced by a simple integration. 

Giauque and Overstreet* seem to have been the 
first to recognize the need for an approximation 
which would take into account the second order 
“stretching” and “‘interaction”’ terms and would 
at the same time avoid the summation over all 
values of J. In the case of chlorine (a typical '}> 
molecule) they deduced an expression which 
replaced by an integral the summation over the 
rotational levels in the state sum for a given 
vibrational state, and showed that the approxi- 
mation was permissible whenever term-by-term 
summation was not feasible (see Eq. (4) below). 
Their method, however, becomes rather laborious 
if many vibrational states must be taken into 
account—either because the temperature is high 
or because the fundamental frequency of the 
molecule is small; moreover, if in addition to the 
free energy, the entropy and heat capacity are 
desired, not only must the state sum be computed 
but also its first and second derivatives with 
respect to temperature, and the calculation may 
then become very laborious. 

In the present paper we describe methods 
which are applicable not only to 1}> but to *>> 
and to *II molecules, and which in many calcu- 
lations eliminate the need for summation entirely; 
even for cases where these methods are not 
immediately applicable, approximations are de- 
veloped which materially reduce the labor of 


3’ Giauque and Overstreet, J. Am. Chem. Soc. 54, 1731 
(1932). 
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computation. One great advantage possessed by 
the procedure outlined below is that it permits 
an immediate estimate to be made of the un- 
certainty resulting in the calculated thermo- 
dynamic quantity from a given uncertainty in 
any of the spectroscopic constants, while a 
second desirable feature is tha ease with which 
previously calculated numerical values can be 
revised in the light of new spectroscopic data. 
The methods are illustrated by calculations for 
various molecules, and since these examples are 
intended to show primarily the precision of the 
approximations, they are sometimes carried to a 
~ greater number of significant figures than is 
justified by the accuracy of the experimental 
data. With regard to the permissible error intro- 
duced by an approximation, it would seem that 
} an accuracy of 0.0005 in the logarithm of the 
| state sum (corresponding to an accuracy of one- 
twentieth of one percent in an equilibrium 
constant), of 0.001R in a molar entropy and of 
0.005R in a molar heat capacity, would leave 

little to be desired by the physical chemist. 

It may be noted in passing that while the 
discussion has been confined to diatomic mole- 
cules, the methods can be readily extended to 
corresponding calculations for the polyatomic 
case, and that the quantities tabulated at the end 
of this paper will be of use in such calculations 
whenever the spectroscopic data become avail- 

able. 


CALCULATION OF THE FREE ENERGY 
Giauque‘ has shown that if F is the free energy 
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in calories at temperature TJ and pressure P 
atmospheres and £,° is the additive energy 
constant at absolute zero, then 


~(F—E?)/RT =3/2-In M+5/2-In T—In P 
+1n R’+1n Q+Eo/RT, (1) 


where JM is the molecular weight, k is the Boltz- 
mann constant, N is the Avogadro number, h is 
Planck’s constant, R is the gas constant in 
calories per degree, R’ is the gas constant in cc- 
atmos. per degree, Q is the state sum defined by 


(2) 


and E, is the energy in the lowest observable 
rotational-vibrational state. In Eq. (2), >, 
indicates summation over all rotational, vibra- 
tional and electronic levels which make appreci- 
able contributions to the state sum, so that for 
moderate and high temperatures, at least a 
double summation is involved. For a hetero- 
nuclear '>> molecule whose constituent atoms 
have no nuclear spin, p,=2J+1 where J is the 
rotational quantum number and e, in cm™! is 
given by 


én = (V+3)we— (V+ +3)? 


where B,=B)(1—av+-yv’), etc., and v is the 
vibrational quantum number. 

Giauque and Overstreet® have shown that for 
small values of B,hc/kT the summation with 
regard to J for any given v may be replaced by an 
integration: 


v 


=etohe/kT, Dem kT. ee -), (4b) 


where in cm~!= (v+1/2)w.—(v+1/2)?w.x, =kT/B,he, d,=(—2D,/B,)qv, It 
should be noted that Qe is a function of v since q,, d,, and f, are all functions of v. Indicating the 
ground vibrational state by the subscript zero, and setting B, = By/(1+8,v+.v?), the expression on 
the right of (4) may be replaced to a good approximation by 


where 


4 Giauque, J. Am. Chem. Soc. 52, 4808 (1930). 
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The errors introduced by taking d and f for the 
ground vibrational state are usually negligible 
but a correction can be made if necessary as 
explained below. Hence to a sufficiently good 
approximation (since 1/12qo, do, fo, 810, etc., are all 
small compared with unity) 


In 0+ Eo/RT = (In +eohc/kT) +1/4¢q0 
go+Bid+ + (6) 


Of the quantities on the right of (6), go, 81, Be, 
d) and fo are known at once, while (In Qy 
+ehc/kT), 6 and v? are functions only of hew,/kT 
and of x; hence double-entry tables of these three 
quantities as functions of w=hcw,./kT and x lead 
at once to an approximate value of (In Q 
+E,/RT). 

It should be noted that if the molecule is 
homonuclear, a symmetry correction — In 2 must 
be added to the right of Eqs. (4) and (6), while 
for each atom possessing a nuclear spin of j, units, 
a spin correction In (2j,+1) must be added to the 
right of the same equations.°* 


CALCULATION OF THE ENTROPY AND HEAT 
CAPACITY 


The entropy of a gas for temperature 7 and 
pressure P atmospheres is computed as the sum 
of the translational entropy Sz and of the 


° These corrections are only approximate, but are valid 
whenever Bhc/kT is so small that Eqs. (4) and (6) are valid 
for the corresponding heteronuclear nonspinning case. As 
an example, consider a rigid homonuclear !2 molecule with 
j.=}, so that the symmetry and spin correction to the 
logarithm of the state sum is — In 2+2 In 2. Assume that 
the weights (as in hydrogen) are (2J+1) for J even and 
3(2J+1) for J odd. The state sum is then 


(A) 
=0,2,--- 


Mulholland (Camb. Phil. Soc. Proc. 24, 280, 1928), has 
shown that 


while from the known properties of the theta-functions 


Hence (A) may be replaced by 


(B) 


For the term 7/480g? to be <0.0001, g must be >12.5, a 
requirement which makes the last term in (B) negligible. 


rotational-vibrational entropy Spy, where 
S7/R=3/2-In M+5/2-In T7—In P 
+1n R’, (7) 
Srev/R=|n In Q/dT, (8) 


and Q is defined in Eq. (4). From Eqs. (8) and 
(6), the approximate value of the rotational 
vibrational entropy is given by 


Srev/R=Sy/R+1+1n 
+15de?/2+3fo, (9) 
Sy/R=In In Qy/aT, 
so=v?+Tdv?/dT. 


where 


Where necessary, the same symmetry and spin 
corrections must be added to the right of (8) 
and (9) as to right of (4) and (6). 
The heat capacity is given by Cp=5R/2+ Cry 
where 
Crv =TdSprv/dT (10) 


(see Eqs. (4) and (8)), while from Eqs. (9) and 
(10) the approximate value of the rotational 
vibrational heat capacity is given by 


+15d?+6fo (11) 


c=Tds2/dT. 


All of the nine quantities (In Qy+ehc/kT), 
Sy/R, Cy/R, i, v?, 1, 52, C1, Ce have been tabulated 
as functions of x and w=hcw,/kT in Tables I-XI 
for x=0.02 and w=0.6. For (In Qy+ehc/kT), 
Sy/R, and Cy/R the value in each case must be 
obtained from two tables; first the value of the 
quantity in question for the given w and for 
x=0 must be found from Table I, and then the 
difference (for the given w) for x=0 and x=x 
must be found from Tables II, III or IV as the 
case may be; this arrangement was adopted to 
facilitate interpolation. In the case of Tables 
I-VIII the interval is sufficiently fine to permit 
interpolation (involving not more than second 
differences) to yield results exact to a unit or so 
in the last place as printed; in the case of Tables 
IX-XI, all that is needed is a rough graphical 
interpolation. Tables II—-IV only extend to w=6; 
the reason for this is that for w>6 an anharmonic 


where 
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TABLE I, w=hcw./kT, x=0. 


w (In Qy +eohc/kT) Sv/R, Cy/R w (In Qy +eohc/kT) Sy/R Cy/R 
0.55 0.8603 1.6103} | 0.9752 2.6 © 0.0772 0.285 0.5859 
60 .7959 1.525 2.7 0696 .264 5631 
65 .7382 1.448 9655 2.8 .0627 2 
.70 1.376 .9602 2.9 0566 
75 6394 1.3108 9544 3.0 0511 .2083 4963 
.80 1.2494 1923 A747 
85 5576 1.1921 9419 3.2 0416 .1776 4536 
90 5218 3.3 .0376 .1640 4330 
0.95 4890 9281 3.4 .0339 1513 4129 
1.00 4587 3.5 .\ 0307 1396 
1.05 4307 .9130 3.6 0277 .1288 3743 
1.10 4048 .9050 3.7 .0250 1188 13558 
1.20 3584 .8882 3.9 .0205 
1.30 .8703 427 0791 .2796 
1.35 .3001 44 0124 .0670 .2436 
1.40 .2832 8515 4.6 0101 .0568 .2170 
1.45 .2673 8418 4.8 .0083 0481 
1.50 S .2525 8319» 59 \- 0068 0407 .1708 
1.55 .2386 8217 5.2 0055 0344 .1508 
1.60 .2255 8114 5.4 .0045 0290 
1.65 .2132 .8010 5.6 .0037 0245 .1168 
1.70 .2017 .7904 5.8 .0030 .0207 1025 
1.75 .1909 .7796 6.0 .0025 0174 .0897 
1.80 .1807 .7687 6.4 .0016 0123 0682 
1.85 1711 .7577 6.8 0011 , 0087 0516 
1.90 .1620 .7466 72 .0061 .0387 
1.95 7.6 .0005 .0043 0.289 
2.00 1454 7241 8.0) 0003 0030 0215 
2.1 .1306 8.4 .0002 .0021 0159 
2.2 1174 .6783 8.8 .0002 .0015 0117 
2.3 .1056 6552 9.2 .0010 
2.4 0951 .6320 9.6 .0007 0062 
2.5 085 .6089 10.0 0045 


oscillator defined by w, and x can be treated 
without sensible error as a harmonic oscillator of 
frequency w,(1—2x); for the same reason Tables 
V-VII only extend to w=4, for w>4 Table VIII 
being used. In all cases entries in the tables were 
computed by term-by term summation of the 
series involved up to the series limit (e.g., v= 24 
for x=0.02). For values of w<0.6, interpolation 
becomes in general impracticable, but for any 
isolated calculations in this range, special ap- 
proximations are available (see Example 4 
below). 

The general method of using the tables can 
perhaps be best illustrated by a numerical ex- 
ample. Assume a heteronuclear '>- molecule with 
nonspinning nuclei, whose energy levels are given 
by Eq. (3) and for which w=hcw,/kT =1.0, 
x=0.02, g=kT/Bohc = 3000, B, = By(1 —0.015v), 


D,=-—10~*B,. Then from Tables I and II, 
(In Qy +eohc/kT) =0.4587 +0.0422; In go = 8.0064, 
1/4qgo=0.0001, B,=0.015, Be=2.25X10-4, d 
=0.0060. From Table V 5=0.687, and from 
Table IX, v?=1.79. Hence 


In Q+ Eo/RT =0.4587 +0.0422 + 0.0001 8.0064 
+0.0103+0.0004+ 0.0060+ 0.0001 = 8.5242. 
The exact value found from Eg. (4) by sum- 
mation up to the series limit v= 24 is 8.5242. To 
compute the entropy, Sy/R=1.0407+0.0975 
(Tables I and III); s; = 1.861 (Table VI); s2=6.96 
(Table X); hence 
Srv/R=1.0407+0.0975+1+8.0064+ 0.0279 
+0.0016+0.0120+0.0003 = 10.1864. 


The exact value (Eqs. (4) and (8)) is 10.1865. To 
compute the heat capacity, Cy/R=0.920/ 
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TABLE IE. (ln Qv+eohc/kT) w, 9 as a 


TABLE III. (S,/R)w,2—(Sv/R)w, 0 as a function of w 


function of w=hcw,./kT and x. =hcw./kT and x. 

w x=0.005 0.010 0.015 0.020 0.025 w x=0.005 0.010 0.015 0.020 0.025 

0.6 0.0169 0.0358 0.0576 0.0837 0.1133 | 0.6 0.0357 0.0778 0.1317 0.2000 0.2695 
0.7 .0143 .0299 .0474 .0677 .0913 | 0.7 .0302 -0649 .1069 .1598 .2206 
0.8 .0123 .0255 .0401 .0566 .0755 | 0.8 .0263 .0557 .0899 -1318 .1822 
0.9 .0107 .0222 .0347 .0485 .0640 | 0.9 .0232 .0488 .0777 .1124 .1534 
1.0 .0095 .0196 -0304 .0422 -0553 | 1.0 .0207 .0434 .0685 .0975 -1318 
1.1 .0084 .0174 .0269 .0372 .0485 | 1.1 .0187 .0390 .0613 0863 ~—-—-«.1153 
12 .0076 .0156 .0241 -0332 .0430 | 1.2 .0171 .0355 .0554 .0775 .1026 
1.3 .0069 .0141 .0217 .0298 0385 | 1.3 .0157 .0325 -0505 .0703 .0923 
1.4 .0062 .0127 .0196 .0269 -0346 | 1.4 .0145 -0299 .0464 .0643 -0840 
1.5 .0057 .0116 .0178 .0244 .0314 | 1.5 .0135 .0277 .0428 -0592 .0770 
1.6 .0052 .0106 .0163 .0222 .0285 | 1.6 .0125 .0258 .0397 .0548 .0710 
id .0048 .0097 .0149 .0203 .0260 | 1.7 0117 0241 .0370 .0509 .0658 
1.8 .0044 -0089 0136 -0186 .0238 | 1.8 .0110 .0225 .0347 .0475 .0612 
1.9 .0040 .0082 .0219 | 1.9 .0103 .0211 .0325 .0445 .0572 
2.0 .0037 .0075 .0115 .0157 0201 | 2.0 .0097 .0199 .0305 .0417 .0536 
in -0031 .0064 .0098 .0134 .0171 | 2.2 .0086 .0177 .0270 .0369 .0473 
2.4 .0027 .0055 .0084 .0114 0145 | 2.4 .0077 .0158 .0241 -0329 .0420 
2.6 .0023 .0047 .0072 .0098 .0124 | 2.6 .0069 .0141 .0216 .0294 .0376 
2.8 .0020 .0040 -0061 .0084 .0107 | 2.8 .0062 .0127 .0194 .0263 .0336 
3.0 0017 .0034 0053 .0072 .0092 | 3.0 .0056 0114 .0174 .0236 .0301 
3.2 .0014 .0030 .0045 .0062 .0079 | 3.2 -0050 .0102 -0156 .0212 .0270 
3.4 .0012 .0025 .0039 .0053 .0068 | 3.4 .0045 .0092 .0140 .0190 .0242 
3.6 0011 .0022 .0033 .0045 .0058 | 3.6 -0040 .0082 .0125 .0170 .0217 
3.8 .0009 .0019 .0028 .0039 .0050 | 3.8 -0036 .0073 -0112 -0152 .0194 
4.0 .0008 -0016 .0024 .0033 -0043 | 4.0 .0032 .0066 .0100 .0136 .0174 
4.2 .0007 .0014 .0021 .0029 .0037 | 4.2 .0029 .0059 .0090 .0122 .0155 
44 .0006 .0012 .0018 0024 .0031 | 4.4 .0025 .0052 -0080 .0108 .0138 
4.6 .0005 .0010 .0015 .0021 .0027 | 4.6 .0022 -0046 .0071 -0096 .0123 
4.8 .0004 -0009 .0013 .0018 .0023 | 4.8 .0020 .0041 .0063 .0085 .0109 
5.0 .0004 .0007 0011 .0015 .0020 | 5.0 .0018 .0036 .0055 .0076 .0097 
5.2 -0003 -0006 .0010 .0013 0017 | 5.2 .0016 .0032 .0049 .0067 .0086 
5.4 .0003 .0005 -0008 0011 .0014 | 5.4 .0014 .0028 .0043 .0059 .0076 
5.6 -0002 -0004 .0007 -0009 0012 | 5.6 .0012 .0025 .0038 .0052 .0067 
5.8 -0002 -0004 -0006 .0008 .0010 | 5.8 0011 .0022 -0033 -0046 .0059 
6.0 0.0002 0.0003 0.0005 0.0007 0.0009 | 6.0 0.0009 0.0019 0.0029 0.0040 0.0052 


+0.1263 (Tables I and IV); c,=2.935 (Table 
VII); c2=21.0 (Table XI); hence 


Cy/R=0.9207 +0.1263+ 1+0.0440+ 0.0047 
+0.0120+0.0005 = 2.108. 


The exact value Eqs. (4), (8) and (10) is 2.109. 

In general, the approximations grow worse the 
smaller w, the larger x, the larger 6 and the 
larger d. For large enough values of d, the error 
introduced by neglecting its variation with v 
may become appreciable; a partial correction 
may be made for this, however, by setting 
d,=d9(1+é6v) and then replacing the dy of (6), 
(9) and (11) by d=d)(1+60). For example, if the 
Constants are those just given except that D, 
=-3X10-B,, then do =0.0180 and d=0,0180(1 


+0.015 X 0.687) =0.01818; if d be used in Eqs. 
(6), (9) and (11), the approximate values are 
8.5371, 10.2129 and 2.137 as compared with the 
exact values 8.5371, 10.2131 and 2.139, while if dy 
is used in the same equations, the approximate 
values are 8.5369, 10.2125 and 2.137 (see 
Examples 2 and 4). 


Example 1*: Cl;,;—(Cl;; at 1000°K 


The constants® are w,.=564.9, w.x«=4.0, B, 
= 0.24295 —0.00170v, D,=—1.83010-7—3.1 


* In none of the examples has the nuclear spin correction 
been added, in Example 1 this would amount to 2 In 6, and 
in example 2, to In 12. 

* Elliott, Proc. Roy. Soc. A127, 638 (1930) as quoted by 
Giauque and Overstreet, reference 3, 
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TABLE IV. (C./R)w, z—(C./R)w, 0 as a function of w 
=hcw./kT and x. 


x=0.005 0.010 0.015 0.020 0.025 


0.0384 0.0936 0.1852 0.2868 0.3213 
.0752 .1416 .2331 .3058 
.0631 .1128 .1863 .2658 
.0544 .0940 1515 .2235 
.0480 .0809 .1263 .1871 

.0430 .0713 .1083 
.0391 .0639 .0950 .1364 
.0358 .0580 .0850 
.0332 .0533 .0772 .1068 
.0310 -0494 .0708 .0969 


0462 d 0888 
0276 0435 
0262 0412 d 0770 
0251 
0241 0376 0691 


0224 d 
0211 0326 d . 
.0309 0547 
.0192 d 0517 
0282 0492 


0177 0271 .0470 
0260 
.0249 
.0156 0238 0411 
0149 0227 0392 


0142 0216 0372 
0135 0205 0353 
0127 0194 
0120 0182 0314 
0112 0171 0294 


0104 0159 0275 

0097 d 0256 

d 0137 0237 
0041 .0083 0127 . 0219 
0.0037 0.0076 0.0116 0.0202 


1 
1 
1 
1 
1 
1 
1 
1 
1 
2 
2 
2 
2. 
2 
3. 
3 
3. 
3. 
3 
4. 
4. 
4. 
4. 
4. 
5. 
5 
5 
5 
5 
6 


F=—9.4X10-'4. Hence,? w=0.80905, 
x=0.00708, g=2873.9, 1/4q=0.00009, 8B, 
= 0.0070, Be =4.9 X 10-5, dp = 0.00433, fo = 0.00002. 
Then (In Q,+ehc/kT) =0.5893+0.0174 (Tables 
I and II); 5=0.850 (Table V); v?=2.3, (Table 
IX). Since the molecule is homonuclear, the 
symmetry correction is — In 2= —0.69315; hence 
the approximate value of (In Q+)/RT) by Eq. 
(6) is 7.8875 while the exact value (Eq. (4)) is 
7.8876. Similarly, Sy/R=1.2387+0.0376 (Tables 
I and IIT); s:=2.121 (Table VI); se=8.1; (Table 
X). Hence, including the symmetry correction, 
Srv/R=9.5707 (exact value by Eq. (8), 9.5709). 
From Tables I and IV, Cy/R=0.9472+0.0393, 
¢,=2.78 (Table VII) ; co =19., (Table XI). Hence 


7In this as in the other examples hc/k was taken to be 
1.43221. 


TABLE V. 0 as a function of w=hcw./kT and x. 


0.005 0.01 0.015 


1.353 1.461 
1.152 

0.936 

654 


557 
479 
415 
362 
317 


.279 


— 


WNNNN 


Crv/R by Eq. (11) is 2.016 while the exact value 
by Eq. (10) is 2.017. 


Example 2: HCl;; at 3000°K 

The constants? are w, = 2994.33, = 53.58, B, 
= 10.4469 —0.31617+0.01195v?, D, = —0.000527 
+0.000012v, F=1.46 10-8. Hence w= 1.42950, 
x=0.01789, go = 200.50, Bi =0.03026, Be 
= —0.000228, d =0.02023(1+0.03775) = 0.02050, 
fo= —0.00034. From the tables (In Qy+eohc /RT) 
=0.2736+0.0231, Sy/R=0.7236+0.0553, Cy/R 
=0.8458+0.0649, 5=0.355, v?=0.64, 51 = 1.036, 
$2= 2.49, €:= 1.670, C2=6.5. Eqs. (6), (9) and (11) 
give 5.6310, 7.1536, and 2.005 as the values of 
(In Q+Eo/RT), Srv/Rand Crv/R, respectively, 
while the exact values (Eqs. (4), (8) and (10)) are 
5.6310, 7.1535 and 2.005. 


Example 3: CO at 3000°K ; 

The constants? are w,= 2167.4, w.«=12.70, 
= 1.843—0.020v, D, = —5.453 X 10-*—6.9 X 
Hence, w=1.0347, x=0.00586, go=1136.5, 6: 


8’ Kemble, J. Opt. Soc. Am. 12, 1 (1926). 

* Snow and Rideal, Proc. Roy. Soc. A125, 462 (1929), as 
quoted by Clayton and Giauque, J. Am. Chem. Soc. 54, 
2610 (1932). 


SOD 


|| 


0.6 1.216 1.279 1o1 
7 0.986 ‘1.031 1.248 
0.9 ‘685 ‘711 ‘821 0 
1.0 ‘582 603 637 
= 499 516 536 583 
431 "445 461 499 | 
375 386 400 432 
296 306 328 | 
253 261 269 280 4 
0139 224 "231 238 246 255 1. 
0132 198 204 211 218 "25 1. 
0126 176 181 187 200 1. 
‘0121 157 162 172 178 2. 
125 128 133 137 142 2. 
0110 100 103 ‘107 ‘110 ‘114 2. 
(0104 ‘080 ‘083 ‘086 ‘089 ‘092 2, 
0098 065 067 069 ‘072 ‘074 2. 
0004 053 ‘054 056 058 ‘060 3. 
043 044 046 048 049 
0087 035 ‘036 037 039 ‘040 
‘ ‘0084 028 029 ‘030 ‘032 033 
0080 023 025 026 .027 
. 0077 0.019 0.020 0.020 0.021 0.022 
0070 
(0066 
\g (0063 
‘0059 
) ‘0055 
i 
tab 


TABLE VI. s, as a function of w=hcw,./kT and x. 
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TABLE VII. c; as a function of w=hcw./kT and x. 


8 


x=0.000 


0.005 


0.010 


0.015 


0.005 


0.010 


0.015 


0.020 


2.834 
2.358 
2.001 
1.724 
1.503 


1.322 
1.171 
1.044 
0.936 


3.023 
2.494 
2.104 
1.804 
1.567 


1.375 
1.216 


3.283 
2.670 
2.231 
1.901 
1.643 


1.436 
1.266 
1.124 


3.706 
2.930 
2.403 
2.024 
1.736 


1.509 


3.734 
3.142 
2.712 
2.384 
2.126 


1.918 
1.745 
1.599 


4.426 
3.577 
3.013 
2.606 
2.297 


2.054 
1.856 
1.692 


5.791 
4.396 
3.526 


6.966 
5.512 
4.355 
3.528 
2.935 


2.515 
2.203 
1.963 
1.771 


(842 : : 1.613 


-625 
569 
519 


433 
363 
306 
.258 
.218 


184 
156 
-132 
112 
0.095 


TABLE VIII. w=hew/kT, x=0. 


5.6 6.0 6.4 
.004 002 002 
025 017 012 
.090 -069 


TABLE X. 52 as a function of w=hcw./kT and x. 


x= 0.000 0.005 0.010 0.015 0.020 


15.70 19.03 25.86 34.43 
12.75 16.10 
9.08 10.88 
6.76 


| 


=0.0109, 62=1.210-4, dy=0.00673. From the | =0.573, v?=1.26, s;=1.508, so=4.40, c,=2.072, 
tables, (In Qy+eohc/kT) =0.4390+0.0107, Sy/R | co=10.s. Eqs. (6), (9) and (11) give 7.4989, 
=1.0092+0.0237, €y/R=0,9154+0.0251, 6 | 9.0994 and 1.979 as the approximate values of 


i 4.280 | 0.6 3.332 1 
3.318 .7 2.854 
99 2.661 8 2.496 
1 2.201 | 0.9 2.217 2.947 
7 CU«datS 1.861 | 1.0 1.992 2.540 
83 || | 1.602 | 1.1 1.808 2.236 
99 1.325 1.397 | 1.2 1.654 1.999 . 
32 1.081 1.173 1.230 | 1.3 1.522 - 1.807 — ‘ 
75 
28 
89 4 : 
55 1 
25 1 
00 1 
78 2 
4? 2 446 459 A74 490 | 2.2 .847 .874 0.902 0.934 0.970 
14 2 A411 | 2.4 .758 .780 .804 831 861 
92 2 315 334 345 | 2.6 -680 .699 -719 .742 -766 
74 2 .266 .274 .283 .292 | 2.8 ~ 611 -627 .645 .664 .685 
60 3 .231 .239 .247 | 3.0 .548 563 .578 595 .613 
49 3 .190 .196 .203 .210 | 3.2 492 505 519 534 550 
40 3 .161 .166 .172 178 | 3.4 441 453 466 478 493 ; 
33 3 -137 141 -146 3S 406 A417 429 442 
27 3 .116 .120 .124 129 | 3.8 .373 384 .396 
22 i 0.098 0.102 0.106 0. 1 10 | 4.0 0.316 0.324 0.334 0.343 0.354 
ue . 
w 4.4 4.8 5.2 6.8 7.6 8.0 
B esi 0.012 -008 -006 .001 .001 .001 .000 
St 0.068 -048 .035 .009 .006 .004 .003 
27 0.250 -196 -153 -052 .039 .029 .022 
0, 
bs 
TaBLe IX. as a function of w=hcw./kT and x. 
‘R 
6, 
13.66 
1) 9.71 
of 7.20 
y, 5.51 
a 4.32 4.71 5.20 5.87 6.96 
re 2.82 3.03 3.29 3.62 4.07 
1.94 2.08 2.23 2.42 2.65 
1.40 1.49 1.58 1.69 1.83 
1.04 1.10 1.16 1.24 1.33 
0.80 0.83 0.87 0.92 0.99 
B 49 | 54 .56 
32 33 35 38 
%. 0.21 0.22 0.23 0.24 0.25 
Br 
as 
54, 
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TABLE XI. ¢2 as a function of w=hcw./kT and x. 


w x= 0,000 0.005 0.010 0.015 0.020 
0.6 30.0 36.7 50.6 80.8 101.3 
0.7 21.6 25.6 32.6 48.4 70.0 
0.8 16.3 18.9 22.8 31.1 44.8 
0.9 12.6 14.4 16.9 21.5 30.2 
1.0 10.0 11.3 13.0 15.8 21.0 
1.2 €.7 7.4 8.3 95 11.6 
1.4 4.7 5.2 5.7 6.3 7.3 
1.6 3.5 3.8 4.1 4.5 5.0 
1.8 2.7 2.9 3.1 3.3 3.7 
2.0 2.1 2.2 24 2.5 2.8 
2.4 1.3 1.4 1.5 1.6 1.7 
2.8 0.9 1.0 1.0 1.1 1.1 
3.2 0.7 0.7 0.7 0.8 0.8 


(In Srv/R and Cry/R, respectively, 
while the exact values (Eqs. (4), (8) and (10)) are 
7.4991, 9.0997 and 1.980. 


Example 4: Brz—Brs; at 1500°K 

The constants” are w,=323.2, w.x«=1.07, B, 
= 0.08077 —0.000275v, D= —2.0310-°. Hence 
w=0.3086, x«=0.00331, q=12966.8, 
= 0.003405, 
=0.00651(1+0.006815). With this value of w, 
the tables are of course inapplicable; but the 
series necessary to evaluate the quantities in- 
volved in Eqs. (6), (9) and (11) can be computed 
without inordinate labor by means of the integral 
approximations listed below. 

The general method is to compute the first few 
terms of the series in question and then replace 
the remaining terms by an integral; the evalua- 
tion of exp [ehc/kT]-Qy may serve as an 
example. The values of exp [(e—e,)hc/kT | 
= exp [—v(1—x)w+v*wx ] for v=1, 2, ---10 are 
computed and added, the result being 3.69403. 
The remaining terms from 11 up are summed by 
the well-known relation 


M 


+(1/12) —As’ J+ (1/24) ] 
—Aw’” 
(12) 


A. R. GORDON AND C. 


BARNES 


where 


Ar’ =f(R+1)—f(R), Ar” =Aryi’—Ar’, 


=Ary" —Ar”’, =Ary” —A,’” etc. 


and N is the series limit. Here 


N N 
f f(y) exp [—y(1—x)w+y’wx ]-dy 
M 


=(1/r)- {exp [—2Mrs+ M?r*]- F(s — Mr) 
—exp [—2Nrs+N?r]- F(s—Nr)}, (13) 


where 
2rs=(1—x)w, 


Fla) =exp [—e?]- f exp [y*]-dy, 


a quantity already tabulated by Lash Miller and 
Gordon.'! Thus by computing f(v) for v= 1, 
M+1, M+2, M+3 and for v=N-—3, N-2, 
N— 1, N (which is sufficient to evaluate the 
difference terms in Eq. (12)), it is possible to 
obtain a numerical value for the sum of the 
remaining terms in the series. 

In the case of bromine, N~150, and for 
1500°K, terms depending on WN are negligible; 
r=0.0319604, s=4.81187; for M=11, s—Ar 
=4.46031, F(4.46031) =0.115163, 2Mrs— 
= 3.25980. Hence the expression on the right of 
(13) =0.13835 and the sum of the terms from 
v=11 to the series limit is 0.15847; therefore 
exp [eke /kT ]- Ov = 3.69403 + 0.15847 = 3.8525; 
the value found by term-by-term summation up 
to v=50 is 3.85258. 

The series involved in Sy/R, Cy/R, 
¢,, and C2 can be evaluated similarly by means of 
Eq. (12), the requisite integrals being as follows 
Writing / for the integral from M to N, E, Eu 
and Ey for 2rsy—r?y?, 2Mrs—M?r?, and 2Nrs 
— N?*r?, and [E], [Em], and [Ey ] for exp [—£] 
exp [— Ey ], and exp [— Ey ], 


SE])-dy=(1/r): F(s— Mr) F(s—Nr)}j, (14) 


© Brown, Phys. Rev. 39, 777 (1932). 


I= ((s— Mr) 


"Lash Miller and Gordon, J. Phys. Chem. 35, 2878-2882 (1931). 
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fF? -(E]-dy = —s?-1,+(s?+ 9) {(s — Mr)- Eu: [Ew 


THERMODYNAMIC QUANTITIES FROM SPECTROSCOPIC DATA 305 


Sy: 


Sy [E]-dy = (28/1) 
S'y*- E-[E]-dy = 


0 


With the aid of the recursion formulae listed 
above, all the integrals necessary for the calcu- 
Jation can be readily computed once J; has been 
found; to preserve numerical accuracy in the 
higher integrals, it is advisable to compute J, to 
two or three decimal places more than are needed 
for the calculation of exp [eohc/kT]-Qry. 

These approximations give (In Qy+ehc/kT) 
= 1.3487, Sy /R=2.2267,Cy/R= 1.0438, = 2.923, 
= 20.5, $:= 6.48, s2=68, c,=7.58, co=1.6X 102; 
d=0.00664, In go=9.4701. Therefore, (In Q 
+E)/RT) = 10.8358, Sev /R= 12.7333 and Crv/R 
= 2.085; the exact values (Eqs. (4), (8) and (10)) 
are 10.8358, 12.7335 and 2.086. The surprisingly 
good agreement in this case in spite of the small 
value of w is due to the fact that x, 6 and d are 
also small. 


Example S: at 3000°K 


The normal electronic state of the oxygen 
molecule is *}° so that the rotational levels fall 
into closely spaced triplets designated by 
Mulliken'? F,, and F3. The weight factors are 
ps =2J+1=2K+3,2K+1, 2K—1 for Fi, Fo and 
F; respectively where K designates the number of 
units of rotational angular momentum. Let e+ 
denote F\(K)—F.,(K) and F;3(K)—F,(K); 
both e+ and «~ vary somewhat with K but 
Without serious error may be taken as constant. 
Let the index + denote constants for F; and — 
constants for F;, constants for F, being without 
index, and for convenience let the zero of energy 


" Mulliken, Phys. Rev. 32, 880 (1928). 


be the state v=0, K =1 for Fy. From the known 
properties of the theta-functions . 
exp [—(K+1/2)*B,he/kT 

—(K+1/2)*- Dyhe/kT ] = (1+3d,/8), (15) 
where q, = kT /B,hce and D, = (—2D,/B,)q.. There- 
fore, since the molecule is homonuclear 
e RT Q = e/kT, e9/4% 


+Qnte } ](16) 
where Qp is defined in Eq. (4). Hence to a good 
approximation 

In Q+ Eo/RT = —In 2+(6)+2/qo+In A, (17) 
Srev/R ~—In2+(9)+5a, (18) 
Crv/R = (11)+ca, (19) 


where (6), (9) and (11) indicate the expressions 
on the right of these three equations with the 
values of go, 81, 82, and dy or d for Fs, and A, sy and 
C4 are defined by 


A = (qo+/qo) (a/qot)!) 4-1 
+ (qo~/qo) (1 — 
$s =In In (20) 
The term 2/go(=9/4qo—1/4go) in Eq. (17) arises 
from the fact that the lowest rotational vibra- 
tional state has J=1, i.e., K=1 for F, while the 


term — In 2 in Eqs. (17) and (18) is the sym- 
metry correction for a homonuclear molecule. 
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The constants" are w.= 1576.78, w,x« = 11.37, 
B,*t = B,=B,~ = 1.4375(1 — 0.01112), D,*=D, 
=D,~ =3.32X10°B,, «*=—-—2.11, e~=—1.90. 
Hence w=0.75276, x=0.00721, qo=1457.15, Bi 
=0.0111, 1.23 X 10-4, d = 0.00968(1+0.01118) 
=0.00978. From the tables (In Qy+ehc/kT) 
= 0.6369+ 0.0193, Sy/R=1.3073+0.0415, Cy/R 
=0.9541+0.0452, 5=0.947, v?=2.84, s;=2.335, 
S2= 9.77, =3.035, co = 23.¢. In the case of oxygen, 
however, there is another complication due to 
the fact that the first excited electronic level, 
which is !}°, is only 13,123 cm~! above the nor- 
mal level. Hence for this temperature, it is 
necessary to add to the expression on the right of 
(16) one-half the product of the expression on the 
right of Eq. (4) into exp [—Aohc/kT] where 
Ay=13,123, the rotational and vibrational con- 
stants for (4) being those for the first excited 
level, viz., = 1426.9, w,’x’ =11.9, B’ =1.3912. 
To a first approximation, the contribution 
from the !} state may be allowed for in Eqs. 
(17), (18) and (19) by adding a term exp 
[ —Aohc/kT ] to the expression for A on the right 
of (20). This of course ignores the change in the 
rotational-vibrational constants in the excited 
state as compared with the normal. If Qy’ denote 
the value of Qy for w= 1426.9hc/kT =0.6812 and 
x =0.00834, this variation can be taken into 
account by adding to the right of (17), (18) and 
(19) respectively 
(21) 
r(Aghc/kT)?- 


A. R. GORDON AND C. BARNES 


where r= | (Qv’— Qv) /Qv + (Bo— B’)/B’} /3. Since 
(In Qy’+e/hc/kT) for w=0.6812 and x = 0.00834 
is 0.7305 while (In Qy+ehc/kT) is 0.6562, the 
correcting terms (21) are 0.00007, 0.0004 and 
0.0028 respectively. 

Including the contributions from the '} state, 
In A=1.09988, s,=1.10321 and c,=0.0248; 
taking account of the corrections (21), Eqs. (17), 
(18) and (19) give 8.3697, 10.0909 and 2.084 as 
the approximate values of (In Q+ Eo/RT), Srv/R 
and Crv/R, while the exact values (Eq. (16)) are 
8.3699, 10.0911 and 2.085 respectively.'4 


Example 6: NO at 3000°K 


In a *II molecule such as nitric oxide, the 
rotational levels fall into two groups, 7ITI\/2 and 
*TI3,2, the energy levels being given by 


+A(J+1/2)+B,(J+1/2)? (22) 
+C(J+1/2)3+D(J+1/2)4, 


where B,=By(1—av+-yv?) = By/(1+810+ 620") 
and J=1/2,3/2,5/2--- for J=3/2,5/2-:- 
for *II3,. The weight factors (excluding nuclear 
spin weight) are 2(2J+1), the factor 2 arising 
from the A-type doubling of the rotational levels. 
The terms involving A and C as determined by 
Jenkins, Barton and Mulliken" are of much less 
importance than the D term, so that 


b-exp [—e,hc/kT ]~2(2J +1) -exp [—(J+4)*Bohc/kT (23) 
—(J+43)*- Che/kT Dhc/kT 


For summation over all positive integral values of » 


~1/b—1/6, = (1/2) (a/b). (24) 


Hence, if constants for IT, be indicated by +, and those for 7II3,, by —, 


13 Mecke and Wurms, Zeits. f. Physik 61, 37 (1930). 
Mecke and Baumann, Zeits. f. Physik 73, 139 (1931). 

4 By using the values of the constants as given in 
International Critical Tables, and by normalizing to the 
state K=1 for F; (which involves adding —2.13hc/kT to 
the right of Eq. (17)), the approximations for 3000°K give 


59.459, 67.848 and 9.108 as the values of (Eo°—*)/T, S 
and C,° respectively. Johnston and Walker (J. Am. Chem. 
Soc. 55, 187 (1933)) found 59.458, 67.848 and 9.108 by 
summation. 

16 Jenkins, Barton and Mulliken, Phys. Rev. 30, 150 
(1927). 
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where +a, +¢,+d,+3d,?), =(—3C/4B,) 


X (1qv)', dy =(—2D/B,.)-q» and ea is the doublet separation. If symbols without index denote the 
quantities obtained by averaging the corresponding quantities for *IIj/. and *II3,2 then to a good ap- 


proximation 


In 2+(1n Oy +eohc/kT) +1/qot+1n go+Bi 


24+ Sy/R+14+1n Bese+ +3¢/2 + 2dy+ 15do?/2 
+ Boece — do/4+ 30/4 + 2do+ 15de?+c,4, 


5do?/2+1n A, (26) 


(27) 


(28) 


where A = qo*/qo+ (qo~/qo) {1 — (2/qo) } 


sa=In A+T0 In 4/dT, 
= T0s,/0T, 


and the quantities 3, s1, se, c1, and are as 
defined under Eqs. (5), (9) and (11). If necessary, 
the dp of Eqs. (26), (27) and (28) may be replaced 
by d=d)(1+60). 

The constants” are: for w,= 1906.54, 
= 14.424, B,=1.675—0.0178v, A =0.0071, C 
=0.000106, D=—5.06X10*; for 
=1906.43, w.«=14.454, B,=1.724—0.0187v, A 
=0.0070, C=0.000010, D=—8.71X10; 
=124.4. Hence w=0.91016, x=0.007574, qo 
= (1250.53+ 1215.00) /2=1232.76, 8,=(0.01063 
+0.01085)/2=0.01074, 
= —0.00010, Co = — (0.00297 + 0.00027) /2 
= —0.00162; do = (0.00756+0.01228) /2 =0.00992 
while d=0.00992(1+0.02155)=0.01007; In A 
=0.66354, s,=0.69269, c,=0.0009. From the 
tables, (In Qy-+eohc/kT) =0.5149+0.0163, Sy/R 
=1.1279+0.0356, Cy/R=0.9337+0.0385, 
=0.713, v?=1.78, 5;=1.823, s2=6.2, ¢,=2.456, 
(= 15.,. Eq. (26) gives 9.0220 as the value of (In 
Q+£E,/RT), the exact value from Eq. (25) being 
9.0221; similarly, Eqs. (27) and (28) give 10.0750 
and 2.022 as the approximate values of Srv/R 
and Cry/R, while from (25), the exact values! 
are 10.7052 and 2.024. 

© With the values of the universal constants as given in 
International Critical Tables, Eqs. (26), (27) and (28) lead 
to the results 60.566, 68.877 and 8.984 for (E\°— F*)/T, S® 


(29) 


Schmid, Koening and von Farkas” in their 
analysis of the rotational structure of the nitric 
oxide spectrum found B,=1.69—0.0175v for 
B,=1.73—0.0190v for *II32, C= —0.0015, 
€a= 120. If these values be used, higher terms in 
the expansion of exp [—(J+1/2)*-Chc/kT] 
must be included and Qe of Eq. (24) will become 
qv‘ The terms 
in a and d will be omitted from Eqs. (26), (27) 
and (28), and co must be replaced by €=c9(1+-60); 
moreover (16/37 —1/2)é? must be added to right 
of (26) and (32/3r—1)é? must be added to the 
right of (27) and (28). Here go=(1239.44 
+1210.79)/2=1225.12, 8:=(0.01036+0.01098) 
/2=0.01067, 1.14 10-4, Co = (0.04154 
+0.04010) /2=0.04082, &=0.04082(1+0.016%) 
=0.04129, In A=0.66446, s,s=0.69271, Ca 
=0.0009. Hence the approximate values of (In 
Q+E./RT), Srv/R and Crv/R are 9.0514, 
10.7463 and 2.036, while the exact values (Eq. 
(25)) are 9.0516, 10.7470 and 2.040. 


and C,° respectively. Johnston and Chapman (J. Am. 
Chem. Soc. 55, 153 (1933)) find by summation for 3000°K 
60.567, 68.877 and 8.982; for 75°K, the approximations 
give 32.747, 40.029 and 7.806, as compared with their 
32.747, 40.028 and 7.806. 

1 Schmid, Koenig and von Farkas, Zeits. f. Physik 64, 84 
(1930). 
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The free energy of steam is computed from the moments 
of inertia determined by Mecke and Baumann, and from 
the fundamental frequencies of Mecke and of Plyler, for 
temperatures up to 3000°K. From these numbers, the free 
energy of carbon dioxide (computed from Martin and 
Barker’s data), the known free energies of carbon monoxide 
and hydrogen, and the thermal data of Rossini, the 
equilibrium constant for the water-gas reaction is obtained. 
The calculated equilibrium constants agree with the best 
experimental values for the temperature range 700°- 
1259°K within the limits of error of the experiments. The 
equilibrium constant for the dissociation of steam into 
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diatomic oxygen and hydrogen is also tabulated, but an 
unambiguous comparison with experiment is not possible 
for this reaction, since steam may also dissociate into 
hydroxyl and hydrogen and the data necessary for a 
reliable calculation of the hydroxyl dissociation are not 
available. A calculation for the producer-gas equilibrium 
predicts greater dissociation of the CO, into CO and 
graphite than is found by experiment; this discrepancy can 
be accounted for if graphite is assumed to have a ‘‘zero 
point” entropy in excess of that given by the third law of 
thermodynamics. 


N a recent paper! the equilibrium constant for 

the water-gas reaction 
was calculated from spectroscopic data. In that 
calculation, Mecke’s older values? for the mo- 
ments of inertia of the water molecule (0.98, 
1.80, 2.8010-) were used and the resulting 
Rin K were in general about 0.2 calory per 
degree less than the experimental values of 
Neumann and Koehler.’ Since that time, two 
independent determinations of the rotational 
constants have appeared, the first by Lueg and 
Hedfeld* and the second by Mecke and Bau- 
mann,°* and as a result a more accurate calcu- 
lation is now possible for reactions involving 
steam. 

Lueg and Hedfeld give for the moments of 
inertia A =0.97, B=2.13, C=3.07 X10-”, while 
Mecke and Baumann give A =1.09, B=1.907, 
C=2.985 X10-*. In the calculation of the ‘‘free 
energy”’ and entropy for moderate and high 
temperatures, it is the effective moment of 
inertia J = (ABC)! which is the significant quan- 
tity, and either set of rotational constants leads 
to approximately the same result. In the follow- 


1 Gordon and Barnes, J. Phys. Chem. 36, 1143 (1932). 

2 Mecke, Trans. Faraday Soc. 26, 214 (1930). 

3 Neumann and Koehler, Zeits. f. Elektrochemie 34, 218 
(1928). 

4 Lueg and Hedfeld, Zeits. f. Physik 75, 512 (1932). 

5 Mecke and Baumann, Phys. Zeits. 33, 833 (1932). 
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ing calculations, the values used were B= 1.907 
x10-", C=2.985xK10-°, A=C—B=1.078 
X10-, giving J=1.83110-; had Lueg and 
Hedfeld’s moments of inertia been selected, / 
would have been 1.851 and the entropy 
and —(F°—£,°)/T (Table I) would have been 
0.033 greater throughout. 

As far as the vibrational contributions to the 
entropy and free energy are concerned, there is 
still an uncertainty as to the third fundamental 
frequency of the water molecule. Mecke® gives 
6=1600, v(x) =3740, v(c)=3850 while 
Plyler? gives 6=1597, v(x) =3742, =5309 
cm-!, For any temperature less than 1000°K, 
either assignment leads to the same value of 
— (F°— E,’)/T within 0.01 cal. /deg. ; for tempera- 
tures greater than 1000°, however, the contri- 
bution of v(c) becomes appreciable, and for this 
reason —(F°—E,°)/T in the range 1100°-3000° 
has been computed both for Plyler’s and Mecke's 
frequencies, the values based on the latter being 
indicated by the star. 

The free energy is obtained from the familiar 


relation 
—(F°—E,°)T 
= —7.267+5R/2-In 7+3R/2-In 18.016 
+3R/2-In /h?—R-In 2+R:In Ov. 


6 Mecke, Zeits. f. Physik 64, 173 (1930). 
7 Plyler, Phys. Rev. 39, 77 (1932). 


| 
| 
t 
f 
n 
| 
err 
col 
N,: 
the 
tio) 
| of ¢ 
of 


FREE 


ENERGY OF STEAM AND CARBON DIOXIDE 


TABLE I. The free energy of sleam. 


T°’K 298.1 300 
—(FP°—E,°)/T 37.226 37.28 
T°K 800 900 
—(F°— /T 45.19 46.18 
T°K 1400 1500 
50.00 50.62 
50.03* 50.66* 

T°K 2400 2600 
/T 55.08 55.88 


400 500 600 700 
39.57 41.36 42.83 44.09 
47.0 47.8 48.64 49.34 
47.07* 47.89* 48.66* 49.36* 
1600 1800 2000 2200 
51.21 52.30 53.30 54.22 
51.26* 52.37* 53.38* 54.33* 
2800 3000 

56.64 57.35 


56.80* 57.53* 


in which J = 1.831 10-", —R In 2 is the Ehren- 
fest symmetry correction, and Q, stands for the 
product of the state sums for three harmonic 
oscillators corresponding to the three funda- 
mental modes of vibration of the molecule. The 
resulting values® of —(F°—£,°)/T are given in 
Table I. The nuclear spin contribution of the two 
hydrogen atoms (2R In 2=2.754) has not been 
included and must be added to the entries if they 
are to be used in conjunction with Giauque’s 
table for hydrogen.® The entropy follows at once 
from the relation S°= —dF°/d7; for 298.1° (in- 
cluding the vibrational entropy 0.008, and the 
nuclear spin entropy 2.754) it is 47.93 cal. /deg. In 
spite of the uncertainty as to the change in the 
spectroscopic constants in the higher rotational- 
vibrational states, the entries in Table I are 
probably accurate to two or three hundredth’s of 


* The starred entries in Table I are computed from Mecke’s fundamental frequencies, unstarred from those of Plyler. 


a calorie per degree up to 1000°, with a con- 
siderably larger uncertainty at the highest tem- 
peratures. 

Of the other gases involved in the water-gas 
reaction, the free energies of hydrogen® and 
carbon monoxide” are known, and an approxi- 
mate calculation for carbon dioxide can be made 
from the data of Martin and Barker''; they give 
I=70.8 X10-”, 6= 667.5, v(s) =1286, v(a) = 2350 
cm~', the deformation frequency 6 being doubly 
degenerate. There is an additional complication 
due to the interaction of v(s) with 6 which tends 
to remove the degeneracy of some of the 6 levels. 
The error introduced by neglecting this inter- 
action, even at 1500°K, is not serious, and for 
lower temperatures is quite negligible. In com- 
puting the entries of Table II, this interaction 


TABLE II. The free energy of carbon dioxide. 


was therefore ignored, but a rough idea of the 


°K 298.1 300 400 500 600 700 800 

—(P°—E,°)/T 43.575 43.62 45.85 47.69 49.28 50.68 51.94 
°K 900 1000 1100 1200 1300 1400 1500 
—(P—Ef)/T 53.10 54.16 55.15 56.08 56.95 57.77 58.54 


error introduced thereby can be obtained from 
considering the interaction of the levels n;=2, 
n.=0, n.=0 and n;=0, n,=1, n.=0. Ignoring 
the interaction, the contribution to the vibra- 
tional state sum is that of three coincident levels 
of energy 2X 667.5 = 1335 cm~! and a single level 
of energy 1286 cm-!. In reality, there are four 


* In the calculations, the values of the universal constants 
are those given in International Critical Tables, Vol. I. 
*Giauque, J. Am. Chem. Soc. 52, 4816 (1930). 


“Clayton and Giauque, J. Am. Chem. Soc. 54, 2610 
(1932) | 


distinct nondegenerate levels” of energy 1286, 
1316, 1340 and 1388 cm™; as far as the calcula- 
tion of —(F°—E,°)/T for 1500° is concerned, 
however, the use of the actual energy values for 
these two quantum states would only change the 
entry in Table II by 0.003 cal./deg. The values 
of —(F°—E,°)/T, computed in the usual way, 
are given in Table II for temperatures up to 


" Martin and Barker, Phys. Rev. 41, 291 (1932). 


” Badger and Sho-Chow Woo, J. Am. Chem. Soc. 54, 
3523 (1932). 
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1500°K, the corresponding entropy for 298.1° 
being 51.09. The numbers, even for the highest 
temperature of the table, are probably reliable to 
a few units in the second decimal place. 


THE WATER-GAS EQUILIBRIUM 


From Rossini’s values for the heats of com- 
bustion of hydrogen and carbon monoxide" 
(68,313 and 67,623 cals., respectively) and 
Osborne, Stimson and Fiock’s value” for the 
latent heat of vaporization of water (2439 joules 
per gram = 10,504 cal. per mole at 25°C), AHags 
for the reaction CO+H2O0=CO2+H: is —9814 
cal. Since the entropies of hydrogen® and carbon 
monoxide” are 33.98 and 47.31 for 298.1°, 
ASoog.1° = — 10.17, and therefore A Fogs. 1° is — 6782 
cal. From Tables I and II and the values of 
— (F°— E,°)/T at 298.1° for carbon monoxide and 
hydrogen, AEo® = — 9636 cal. The resulting RIn K 
=R pco- px20/pcog: for the range 600° 
1500° are given in Table III; in arriving at these 


TABLE III. The equilibrium constant for the water-gas 
reaction. 


K=bco: bu:0/ Pr:- 


T°K 600 700 800 900 1000 
—-RinK 6.49 4.32 2.74 1.54 0.62 
T°K 1100 1200 1300 1400 1500 
-RinK -O11 -0.70 -1.18 -158 -—1.91 


values, entries in Table I based on Plyler’s 
frequencies were used. 

Neumann and Koehler* found as the mean of a 
large number of measurements at 1259°, K = 1.61, 
i.e., R In K=0.947; the value interpolated from 
Table II for this temperature is 1.00, i.e., 
K=1.65. With four exceptions, their results for 
700° to 1259° give R In K which lie within +0.12 
cal. /deg. of the numbers obtained from Table ITI. 
For the range 1100°-1259° (measurements at 
seven different temperatures) the mean deviation 
of the experimental R In K from the calculated is 
— 0.05; for the range 900°-1100° (measurements 
at eight temperatures) the mean deviation is 

—0.04; and for the range 700°-900° (measure- 


3 Rossini, Bur. Standards J. Research 6, 1 (1931). 
4 Rossini, Bur. Standards J. Research 6, 37 (1931). 
% Osborne, Stimson and Fiock, Bur. Standards J. 


Research 5, 411 (1930). 
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ments at five temperatures) the mean deviation 
is +0.06. 


Emmett and Schultz" obtained values of the 


equilibrium constant for this reaction from their 
data on the reduction of cobalt oxide by carbon 
monoxide and by hydrogen. The values of RIn K 
from their measurements at 723°, 788° and 843° 
are —3.95, —2.90 and —2.14, respectively, as 
compared with —3.91, —2.90 and —2.18 from 
Table III. 

Giauque and Ashley” have shown that the 
spectroscopic data lead to a value for the free 
energy of formation of steam in agreement with 
the best measurements on the mercuric oxide 
equilibria, and also in agreement with some 
calculations of their own for the Deacon reaction. 
It has been shown here that the same data lead to 
an equilibrium constant for the water-gas re- 
action which agrees with the most careful 
measurements of that quantity within the limits 
of error of the experiments. 


THE THERMAL DISSOCIATION OF STEAM 


From Johnston and Walker’s value'® for the 
entropy of oxygen at 298.1°, viz., 49.02, ASoos.1° 
for the reaction He+}302=H2O(g) is — 10.56; 
combining this with Rossini’s value” for AHo9s.1 
(—57,809 cal.), —54,661; therefore, 
from the known free energies of hydrogen® and 
oxygen,'® AEy® = —57,121. The resulting R In K 
=R In (pu20) / (pu) - (poe)! are given in Table IV 
for temperatures!® from 1300° to 3000°. As in the 


TABLE IV. The equilibrium constant for the “oxygen” 
dissociation of steam. 


K = (p20) (bo.)!. 


T°K 1300 1400 1500 1600 1800 = 2000 
RinK 32.33 29.05 26.18 23.67 1948 16.12 
T°K 2200 2400 2600 2800 3000 
RinK 13.36 11.05 9.07 7.38 5.90 


1% Emmett and Schultz, J. Am. Chem. Soc. 52, 1782 
(1930). 

17 Giauque and Ashley, Phys. Rev. 43, 81 (1933). Their 
value of AF for steam is —54,669, while the value obtained 
here is —54,661. The slight discrepancy arises from the 
fact that they used 10,490 as the latent heat of vaporization 
and 47.92 as the entropy. 

18 Johnston and Walker, J. Am. Chem. Soc. 55, 172 
(1933). 

19 The oxygen numbers are subject to some uncertainty 
at the higher temperatures owing to the possible existence 
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case of the water-gas reaction, the numbers in 
Table I used in computing Table IV were those 
based on Plyler’s frequencies. 

A comparison of Table IV with the experi- 
mental data for the dissociation of steam is 
subject to some uncertainty. If the experimental 
results be accepted at their face value, the 
agreement between the R In K obtained by 
various workers” for the temperature range 
1300°-2200°K and the numbers obtained from 
Table IV is surprisingly close—in general within 
+0.3 cal./deg. Such agreement must be to some 
extent fortuitous, however, since Bonhoeffer and 
Reichardt®! have shown, qualitatively at any 
rate, that the ‘“‘hydroxyl”’ equilibrium HxO =OH 

>He is at least as important as the ‘‘oxygen”’ 
equilibrium computed in Table IV. Unfortu- 
nately a reliable calculation of the hydroxyl 
dissociation is not possible at present owing to 
the uncertainty as to the value of AE, for the 
reaction; until this point is settled, a comparison 
of the experimental results with Table IV is 
largely guess-work. 


THE PROpDUCER-GAS EQUILIBRIUM 


From Rossini’s'* heat of combustion of carbon 
monoxide (67,623 cal.) and Roth and Naeser’s 


a'A level; see Lewis and von Elbe (J. Am. Chem. Soc. 55, 
511 (1933)). Such a level makes a contribution to —(F° 
—E,°)/T of approximately 2R/3 exp [—17,200/RT], i.e., 
Johnston and Walker’s values of —(’—E,°)/T will be 
increased by 0.004 at 1500°K and by 0.074 at 3000°K; the 
corresponding R In K will be decreased by one-half these 
amounts. 

Int. Crit. Tab., Vol. VII, p. 231. 

*1 Bonhoeffer and Reichardt, Zeits. f. physik. Chemie 
139, 75 (1928). Bonhoeffer and Reichardt’s calculation of 
Kou=(Pu.0)/(Pu:)':(pou) is in’ error, since they did 
not allow correctly for the multiplicity of the rotational 
levels of the hydroxy! molecule. By using the spectroscopic 
constants given by Weizel (Bandenspektren, Leipzig, p. 
327 (1931)) the values of —(F°—E,°)/T for hydroxyl are 
48.35, 50.44, 52.11 and 53.50 for 1500°, 2000°, 2500° and 
3000° respectively. The calculation is similar to that for 
NO discussed by Gordon and Barnes (J. Chem. Phys. 
1, 297 (1933)) except that the I-doublet is inverted, with 
the result that the lowest observable state has v=, 
J=3/2 for *II,,.; the numbers just given include the item 
R\n 2 which allows for the A-type doubling of the rotational 
levels but do not include the spin contribution of the 
hydrogen atom. If Bonhoeffer and Reichardt’s estimate 
Cal. for the reaction }H:+OH=H,0 be 
accepted, then R In Koy for the four temperatures is 27.2, 
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TABLE V. The equilibrium constant for the producer-gas 
reaction. 


K 


T°K 600 700 800 900 1000 
RinK —-—25.91 -—16.10 —8.77 —3.10 +1.42 
T°K 1100 1200 1300 1400 1500 
RinK +5.09 +8.14 +10.69 +12.87 +414.75 


heat of combustion of graphite” as corrected by 
Clayton and Giauque,” viz., 94,240 cal., AHogs.1 
= 41,006 cal. for the reaction C(8 graphite) + 
=2CO. If the “third law” entropy of graphite 
(1.39 at 298.1°) be accepted, then ASo9s.1° = 42.15 
and AFo3.;° is 28,441 (corresponding to a free 
energy of formation for CO: at this temperature 
of —94,443 cal.) From Clayton and Giauque’s 
tables for carbon monoxide and graphite, and 
Table II, AEy° = 39,352 cal.; the resulting R In K 
=R In (pco)?/(peoe) are given in Table V for 
temperatures from 600° to 1500°K. 

The experimental R In K for the temperature 
range 1073°-1473°, obtained from the measure- 
ments of Rhead and Wheeler,” of Jellinek and 
Diethelm,”* of Dent and Cobb” and of Becker** 
are quite self-consistent but are uniformly about 
0.5 cal. /deg. less than the corresponding numbers 
obtained from Table V. Of the quantities used in 
computing the table, there is no reason to suspect 
the heats of combustion or the spectroscopic 
entropies and free energies of the oxides of carbon. 
The position is different with respect to the 
graphite numbers, since it is possible that there is 
a “zero-point” entropy not included in the 
extrapolated integration to absolute zero. If ad 
hoc, for the sake of argument, a zero point 
entropy of 0.5 cal./deg. be assumed for graphite, 
the entropy at 298.1° will be 1.89, and all the 
entries in Clayton and Giauque’s table of 
—(F°—E,°)/T will be increased by 0.5. AFo9s.1° 
for the producer-gas reaction will then be 28,590 


16.0, 9.3 and 4.8 respectively. Since AEo® is possibly in 
error by several kilo-calories, these R In Koy are uncertain 
to the extent of several calories per degree. 

22Roth and Naeser, Zeits. f. Elektrochemie 31, 461 
(1925). 

23 Rhead and Wheeler, J. Chem. Soc. 99, 1140 (1911). 

24 Jellinek and Diethelm, Zeits f. anorg. Chemie 124, 203 
(1922). 

25 Dent and Cobb, J. Chem. Soc. 1929, p. 1903. 

26 Becker, J. Iron Steel Inst. (London) 121, 337 (1930). 
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cal., all the RIn K in Table V will be decreased by 
0.5 cal./deg., and the free energies of formation 
of the two oxides will be —32,850 cal. and 
— 94,290 cal. for CO and COs, respectively. In a 
recent survey of this and other technical gas 
reactions, Chipman” gives 28,570 as the ‘‘best”’ 
experimental value of AFo9s.:° for the producer- 
gas reaction; the result obtained by assuming a 
zero point entropy for graphite of 0.5 cal. /deg. is 
in very close agreement with his value. 

A natural way to settle the question as to the 


27 Chipman, Ind. Eng. Chem. 24, 1013 (1932). 


proper value of the entropy of graphite would be 
to compute it from the entropy of carbon vapor, 
the vapor pressure of graphite, and the heat of 
sublimation. Unfortunately, the requisite experi- 
mental data?* are not known with sufficient 
accuracy to permit such a calculation. The 
results for the producer-gas reaction suggest, 
however, that in other reactions involving 
graphite the possibility that the ‘third law” 


entropy is too small should be borne in mind. 


28 See a recent paper by Vaughan and Kistiakowsky, 
Phys. Rev. 40, 457 (1932). 
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Cross-Activation in the Unimolecular Decomposition of Mixtures of Gaseous Methyl’ 
and Ethyl Ether 


E. W. R. Steacie, Physical Chemistry Laboratory, McGill University, Montreal 
(Received March 13, 1933) 


On account of the disagreement between previous 
investigations of Kassel and the writer, the decomposition 
of mixtures of methyl and ethyl ether has been rein- 
vestigated with greater accuracy. The cross-activation 


efficiency has been found to be considerably greater than 
that previously reported. It is, however, decidedly smaller 
than that for collisions between molecules of the same kind. 


INTRODUCTION 


N previous communications! it was reported 
that mixtures of gaseous methyl and ethyl 
ether decomposed at rates determined solely by 
the partial pressures of the two ethers, and hence 
the efficiency of energy transfer in cross collisions 
of the second kind was apparently zero. Kassel? 
repeated the work and found efficiencies between 
66 and 153 percent, depending on the pressure. 
He therefore concluded that, while the actual 
values were somewhat uncertain, the cross-activa- 
tion efficiency must be close to 100 percent. 

The present paper describes a repetition of the 
work with considerably greater accuracy, in an 
attempt to settle the question of the discrepancy 
between the two investigations. 


SOURCES OF ERROR 


Variations in temperature are almost invari- 
ably the most important source of error in 
investigations of the kinetics of homogeneous gas 
reactions. This is especially true in the case of 
unimolecular reactions on account of their high 
temperature coefficients. In the case of the ether 
decompositions, the rates vary from 10 to 15 
percent per degree difference in temperature. 

In the previous investigations of Kassel and 
the writer ordinary electric furnaces were em- 
ployed. In both cases the average temperature 


*Steacie, J. Am. Chem. Soc. 54, 1695 (1932); J. Phys. 
Chem. 36, 1562 (1932). 
_ * Kassel, J. Am. Chem. Soc. 54, 3641 (1932). I am 
indebted to Dr. Kassel for communicating his results to me 
Prior to their regular publication. 


was probably constant to within about 0.5°C. 
This, however, is equivalent to about 6 to 7 
percent error in the rates of reaction, and will 
largely mask the effect under investigation. Any 
attempt to obtain greater accuracy must there- 
fore involve considerably better temperature 
control. 


EXPERIMENTAL PROCEDURE 


The reaction velocities were followed in the 
usual way by admitting the reactant or reactants 
to a heated silica bulb connected to a capillary 
manometer, and observing the change in pressure 
at regular intervals of time. Except for the 
temperature bath, the apparatus was identical 
with one which has been used in previous 
investigations.* The necessary temperatures were 
obtained by immersing the reaction bulb in a 
well-stirred bath of molten lead and tin (80 
percent lead). 

The bath consisted of a 12 inch length of 7 inch 
wrought iron pipe, to the bottom of which was 
welded a flat iron plate. It contained about 4 
liters of molten metal. The iron pipe was covered 
with asbestos paper and wound with 16 gauge 
nichrome wire. The bath was packed in a steel 
container with fire-brick and loose asbestos. It 
was stirred vigorously by a wrought iron stirrer, 
the shaft of which passed through a small hole in 
the fire-brick lid. 

Temperatures were measured with a chromel- 
alumel thermocouple in conjunction with a 
Cambridge thermocouple potentiometer which 


3 Steacie, Can. J. Research 6, 265 (1932). 
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could be read to about 0.05°C. The thermocouple 
was checked at the freezing and boiling points of 
water, and at the boiling point of sulphur. 

Rough temperature control was effected by a 
200 cc Pyrex bulb filled with air and immersed in 
the bath, and connected to a thermal regulator of 
the ordinary type. With this the temperature 
could be maintained constant to within about 
2°C. During a run the temperature was controlled 
by hand regulation of rheostats. 

On account of the high heat capacity of the 
contents of the bath, fluctuations in temperature 
were reduced to a minimum, while the high heat 
conductivity of the lead-tin mixture and the 
stirring ensured uniformity of temperature. The 
temperature could be maintained constant to at 
least 0.1°C. No variation in temperature through- 
out the bath could be detected. 

The efficiency of the heat insulation is shown 
by the fact that a temperature of 478°C could be 
maintained by a power consumption of about 350 
watts. On account of the efficient insulation, the 
bath was almost air-tight, and the lead-tin 
mixture did not oxidize to any great extent. 
After about two months continuous operation at 
478°C, only about 5 percent of the lead-tin 
mixture had oxidized. 

Methyl and ethyl ether were prepared and 
stored as previously described. 


EXPERIMENTAL RESULTS 


The relative rates of reaction were inferred, as 
in Kassel’s investigation, by comparing the 
fractional pressure increases at 5 minute intervals 
for 20 minutes in each run. For the first 5 minutes 
readings were taken at 1 minute intervals and 
extrapolated to zero time to give the initial 
pressure. Runs with the separate ethers and with 
the mixture were made in random order. The 
values of the fractional pressure increase in 20 
minutes for various initial pressures are given in 
Tables I, II and III. 

The data are plotted in Figs. 1, 2 and 3. In 
order to obtain lines of small curvature, the 
values of the fractional pressure increase are 
plotted against logio Py rather than against Py 
itself. 

It was considered desirable to draw the curves 
through the experimental points by a compu- 
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TABLE I. Fractional pressure increase in 20 minutes for 
ethyl ether at 478°C. 


Po Po 
cm AP 20/Po cm. AP 20/Po 
35.40 0.2720 15.89 0.2128 
35.34 0.2610 14.40 0.2143 
34.45 0.2669 14.30 0.2080 
33.92 0.2642 13.27 0.1989 
31.49 0.2740 12.98 0.1918 
29.95 0.2530 12.89 0.1850 
28.93 0.2535 10.91 0.1803 
26.65 0.2492 10.20 0.1843 
24.76 0.2355 9.42 0.1632 
23.38 0.2578 9.23 0.1700 
22.90 0.2330 9.06 0.1622 
22.25 0.2375 8.79 0.1801 
20.94 0.2238 7.58 0.1660 
20.53 0.2265 7.55 0.1629 
18.76 0.2242 7.33 0.1498 
18.38 0.2233 7.19 0.1680 
17.74 0.2161 6.88 0.1584 
17.63 0.2318 6.56 0.1645 
17.10 0.2122 6.50 0.1493 
16.42 0.2070 5.70 0.1486 


TABLE II. Fractional pressure increase in 20 minutes for 
methyl ether at 478°C. 


Po Po 

cm AP 20/Po cm AP 2o/Po 
60.23 0.3250 15.46 0.1934 
53.10 0.3130 14.52 0.1665 
52.65 0.3120 13.68 0.1746 
51.80 0.2938 13.37 0.1533 
42.57 0.2800 11.37 0.1530 
35.86 0.2510 10.89 0.1478 


24.61 0.2180 6.95 0.1180 

6.74 0.1317 
21.74 0.2150 6.19 0.1117 
19.07 0.1892 6.09 0.1132 
17.83 0.1926 5.41 0.1182 
17.53 0.1778 4.14 0.1037 
17.49 0.1987 


tational method, so as to avoid any doubts as to 
the objectivity. The curves have therefore been 
drawn by fitting the experimental results to 
empirical equations of the type 


AP»/Po=a+b log Po+c(log Po)’. 


The fitting has been done by the method of 
averages, since the additional labor involved in 
the use of least squares did not seem to be 
warranted by the experimental accuracy. 


35.50 0.2644 9.84 0.1484 

35.13 0.2570 9.54 0.1440 

34.89 0.2455 8.86 0.1420 

ee 30.35 0.2542 8.51 0.1468 

_ 27.10 0.2250 7.82 0.1176 
ay 24.84 0.2136 7.48 0.1257 
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TABLE III. Fractional pressure increase in 20 minutes for a 
one to one mixture of ethyl and methyl ether at 478°C. 


AP 20/Po 


0.2298 
0.2142 
0.2175 
0.2238 
0.2165 
0.2075 
0.1950 
0.2030 
0.1985 
0.2047 
0.1846 
0.1885 
0.1824 
0.1682 
0.1708 


AP 20/Po 


0.1657 
0.1722 
0.1700 
0.1500 
0.1517 
0.1525 
0.1485 
0.1391 
0.1356 
0.1390 
0.1190 
0.1284 
0.1248 
0.1193 
0.1275 


6.70 
6.27 


Leg, PR 


Fic. 2. Methyl ether. 


12 


Log, 


Fic. 3. 1 : 1 mixture. 


The cross-activation efficiencies have been 
calculated by comparing the values of the 
fractional pressure increase in 20 minutes for the 
mixture with those calculated for the separate 
ethers: (A) on the assumption that each ether in 
the mixture will decompose at a rate determined 


by its partial pressure, i.e., on the assumption of 
zero efficiency of energy transfer between unlike 
molecules; and (B) assuming that the rate of 
decomposition of each ether in the mixture will 
correspond to its rate of decomposition when 
alone at a pressure equal to the fofal pressure of 
the mixture, i.e., assuming 100 percent efficiency 
of energy transfer between unlike molecules. 
The results are shown in Table IV. 

Since complications may arise during the later 
stages of the reactions, it is obviously desirable to 
compare the relative rates of reaction at the 
earliest possible stage. Kassel accomplished this 


TABLE IV. Cross-activation efficiencies. 


Fractional pres- 
sure increase 
calculated for 
mixture, 
assuming 


Per- 

cent 

cross 
activa- 


Fractional pressure 
increase in 20 
minutes, observed 100 


percent Zero tion 
Po Ethyl Methyl Mix- effi- effi- 
cm Po ether ether ture 


ciency ciency ciency 
0.125 0.101 0.095 


0.113 
0.137 0.109 0.107 0.123 
0.151 0.119 0.120 0.135 
0.165 0.131 0.133 0.148 
0.180 0.145 0.146 0.163 
0.195 0.161 0.160 0.178 
0.211 0.179 0.175 0.195 
0.227 0.199 0.190 0.213 
0.243 0.221 0.205 0.232 
0.261 0.246 0.222 0.254 
0.279 0.272 0.239 0.276 
0.298 0.300 0.258 0.299 
0.316 0.330 0.277 


ooooo 


S38 


0.323 


| | 
18 p 31.87 14.95 
13 O20 31.78 13.89 
31.65 13.77 
0) ows 25.12 10.76 
24.98 10.72 
13 24.87 9.99 
2 22.07 8.54 
y) 19.51 6.97 
1 17.70 6.85 
0) 17.53 
1) 15.61 
18 15.26 
5 
= 030 
for 
020 
6 | 
3 
4 
ozs 
8 
6 
020 
0 
7 
32 
ols 5 
4.0 
en || 6.3 
7.9 0.113 57 
to 10.0 0.123 58 
12.6 0.135 58 
15.8 0.148 57 
20.0 0.163 54 
25.1 0.178 50 
31.7 0.195 46 
of 39.8 0.213 43 
i 50.2 0.232 40 
in 63.2 0.254 39 
be 
Mean 50 


316 W. 
by determining the fractional pressure increases 
at 5, 10, 15 and 20 minutes, and extrapolating 
back to zero time to get the initial rate of 
reaction. An inspection of Fig. 2 of Kassel’s 
paper, however, shows that there is no uniform 
trend to the results with increasing increments of 
time, and within the experimental error virtually 
the same result is obtained by comparing values 
of AP2/Po, as by comparing the initial rates. 

In the present investigation it was found that 
efficiencies calculated from a comparison of 
fractional pressure increases at 5, 10 and 15 
minutes showed no significant difference from 
those tabulated above. No appreciable error is 
therefore introduced by the use of fractional 
pressure increases at 20 minutes as a criterion of 
the relative rates of reaction. 


DISCUSSION 


As mentioned at the outset, the previous 
investigation of the writer led to the conclusion 
that the cross-activation efficiency in mixtures of 
ethyl and methyl ether was zero, while that of 
Kassel indicated an efficiency of the order of 66 
to 150 percent depending on the pressure. The 
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present investigation gives values between these, 
viz., 40 to 60 percent. It is almost certain that the 
discrepancies between the former investigations 
were due to temperature variations partially 
masking the effect under investigation. The 
values obtained in the present investigation 
should be considerably more reliable, on account 
of greatly improved thermostating. 

In Fig. 3 the dotted curves indicate the values 
of AP»/Po» for the mixture of ethers calculated 
on the basis of zero and 100 percent efficiency of 
energy transfer between unlike molecules. A 
comparison of these two curves with the experi- 
mental points indicates that a fair amount of 
reliance may be placed on the numerical values of 
the cross-activation efficiencies given above. 

It may therefore be concluded: (a) that the 
contention of Kassel is borne out that collisions 
between unlike complex molecules are definitely 
effective in producing activation in unimolecular 
reactions; (b) that there is, however, a definite 
specificity of energy transfer, collisions between 
unlike complex molecules being considerably less 
effective than those between molecules of the 
same kind. 
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The Normal State of the Hydrogen Molecule-Ion* 


B. N. Dickinson, Gates Chemical Laboratory, California Institute of Technology 
(Received March 16, 1933) 


A variation function somewhat similar to that applied by Rosen to H2 is used in the discussion of the normal state of H2*. 


N treating the normal state of the hydrogen 

molecule by the variational method, N. 

Rosen! used an eigenfunction involving the one- 
electron function 


v=P+oy', 


in which y°=(a*/)!e-* is a hydrogen-like 1s 
wave function for nucleus a, with effective 
nuclear charge (in units e) Z=ado, and y’ 
=(a5/r)'r,e-* cos 60, is a hydrogen-like 2p 
function for nuclear charge 2Z, @, being the angle 
between the nuclear axis and the radius vector 
from nucleus a to the electron, and o a parameter. 
In order to see to what extent this method of 
treatment of molecules could be expected to be 
satisfactory in other cases, and whether or not 
pronounced improvement could be obtained by 
varying the ratio of the effective nuclear charges 
in ¥° and y’, we have treated the hydrogen 
molecule-ion (which has previously been well 
investigated) by methods similar to Rosen’s. 

For this purpose we replace Z by Z’ =e«Z in y’. 
Y° and y’ satisfy the equations 


TYp=Z2EwW and 
where 


—h? 


Zé 


822m Ye 


—h? 2Zee* 
Vv? = ’ 
822m Yo 


and H,'= 


and are normalized and mutually orthogonal. 


We write 

as the complete eigenfunction for the molecule- 
ion, with @° and @’ replacing y° and y’ for the 
nucleus b. The Hamiltonian for this problem is 


822m 


The integrals arising from the expression for the 


energy 


are listed below using Rosen’s notation. 


f 


To = 


In= [ 


Ya 


The integrals (Ju, Gor, Jor, and Jio) not given in 
Rosen’s paper were evaluated in elliptic coordi- 
nates and have the following forms in terms of «€ 
and p=aR: 


To, = (8e*/2/(1 — { — 1) 4p? — — 24€} (1 — 8ep?(1 p+ 24€} J, 
RJu = (46!?/(1 — — +1) { (1 — — 4e(1 — p? + + 8e} J, 
= (8€/?/(1 (1 p? +4} { (1 —4ep—4} ], 


RG = — { (1+ + 4(1 +8]. 
*Contribution from Gates Chemical Laboratory, California Institute of Technology, No. 363. 


*N. Rosen, Phys. Rev. 38, 2099 (1931). 
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The energy can be expressed as 


N. DICKINSON 


W=22°E)A /d+2ZEB/d, 


where 


A=1 +Io+oe(1 In +o°e(1 = (2/p)[RFot+ + |, 
B = (2/p) LR Fot+ RGoo— 1 +2RJo+20(RJu— Tu +RGout+RJio) 1) ], 


d = 2[1 +Jn)]. 


This expression for the energy can be mini- 
mized with respect to Z, giving 


= —B/2A, 


which when combined with the expression for W 
gives 
W = — EyB?/2dA. 


RESULTS 


On placing e=1 and minimizing the energy 
relative to Z, o, and p, the values e=1, p=2.50, 
o=0.1605, Z=1.2537, and W=—16.242 volt- 
electrons were obtained. The value given for W 
by Burrau’s? and Hylleraas’ exact treatments 
is —16.30 volt-electrons. Consequently the 
Rosen-type function leads to an error of 0.06 
volt-electron, which is rather unsatisfactory, in 
view of the fact that Guillemin and Zener* ob- 
tained the value W= —16.30 volt-electrons by 
using a two parameter function no more compli- 
cated than the Rosen-type function. 

It was thought that the error of 0.06 volt- 
electron in W might be removed by varying e. On 
carrying out this variation for the values of o, 
Z, and p given above, then varying oa, and finally 


2 An account of this work is given by Linus Pauling, 
Chem. Rev. 5, 173 (1928). 

3 E. A. Hylleraas, Zeits. f. Physik, 71, 739 (1931). 

4V. Guillemin and C. Zener, Proc. Nat. Acad. 15, 314 
(1929). 


varying ¢ again, there were obtained the values 
p=2.50, e= 1.15, o =0.145, 
Z=1.247, Z'=1.434, 1rm=1.06A, 


and W = — 16.257 volt-electrons. The discrepancy 
in the energy has, it is seen, been reduced by only 
one-third by the variation of ¢, being still 0.04 
volt-electron. (The value of 7 obtained is the 
same as that found by Burrau and Hylleraas and 
also found experimentally for excited states of 
He.) We accordingly conclude that the Rosen- 
type function for the hydrogen molecule-ion is 
not so satisfactory as the Guillemin-Zener-type 
function of equivalent complexity, and that no 
pronounced improvement results from a variation 
of « from 1. Change of Rosen’s treatment for the 
hydrogen molecule by introducing « in his func- 
tion and varying it can hence be expected to lead 
to little improvement over his results. 

It must be mentioned that the error in W of 
0.04 volt-electron is only about 8 percent as great 
as that for a function with o =0, i.e., involving 1s 
functions alone, for which Finkelstein and 
Horowitz®> obtained Z=1.228, 7 =1.06A, and 

= — 15.75 volt-electrons, so that very marked 
improvement in the wave function is obtained by 
introducing the Rosen-type term y’. 


5 B. N. Finkelstein and G. E. Horowitz, Zeits. f. Physik 
48, 118 (1928). 
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Lifetimes of Unstable Molecules 


N. RoseEn,* University of Michigan 
(Received February 1, 1933) 


The mechanism of the decomposition of a molecule 
formed by the collision between a molecule and another 
particle is discussed qualitatively. Calculations are then 
carried through to determine the lifetime of a simple linear 
triatomic molecule formed by the collision of an atom witha 


diatomic molecule. The results are applied to the short- 
lived molecule HO: to determine the order of magnitude of 
its lifetime in a number of different states, although it 
does not conform strictly to the idealized conditions for 
which the calculation was made. 


I. 


HERE is evidence that, as intermediate 

steps in certain chemical reactions, mole- 
cules are formed of extremely short lifetime, 
decomposing spontaneously unless they can react 
with other molecules before this happens. In 
some cases the molecules are known to be 
intrinsically stable. The question naturally arises 
as to the mechanism of this unimolecular de- 
composition. While this type of phenomenon has 
been treated by other writers, it is the purpose of 
the present paper to discuss it in slightly greater 
detail than has been done hitherto and from a 
somewhat different point of view. 

Let us first consider the formation of a mole- 
cule from two colliding particles. In order that 
they may form a stable system after the collision 
two conditions must be satisfied. First the force 
between the particles must be such that their 
potential energy as a function of the distance 
between them must have a minimum, that is, its 
graph must have the form of a valley. Then, the 
sum of this potential energy and the energy of 
relative motion of the particles must be less than 
the top of either side of the potential valley. It is 
evident that since the particles were originally 
free and therefore had energy above the outer 
wall of the potential valley there must be 


* National Research Fellow. 

‘For example, Roginsky and Rosenkewitsch, Zeits. f. 
Physik. Chemie B10, 47 (1930), Langer, Phys. Rev. 34, 92 
(1929), Kassel, J. Phys. Chem. 32, 1065 (1928), etc., Rice, 
Phys. Rev. 34, 1451 (1929), etc. The paper of Roginsky and 
Rosenkewitsch treats the problem for the standpoint of the 
Auger effect as is done in the present paper. 


something to take up some energy in order to 
bind the particles to each other. As is well known, 
this can be done by a third particle (triple 
collision), the wall of the containing vessel, etc. 

It must be remembered, however, that it is not 
necessary to remove this energy completely from 
the system. It is sufficient merely to transfer it 
from the degrees of freedom corresponding to the 
relative motion of the two particles to some other 
part of the system. Thus, if the colliding particles 
are both atoms the energy might go to excite the 
electrons in one or both of them, producing a 
molecule in an excited electronic state. If one of 
the particles is a molecule there is an additional 
way of disposing of this energy. It might be taken 
up by the internal degrees of freedom of this 
molecule as vibrational energy. It is this case 
that interests us in the present paper, namely, 
that in which a collision takes place between a 
molecule and an atom or another molecule so 
that some vibrations are excited in the former and 
a binding to the latter thus produced. 

With this picture it is easy to see why the 
molecule, if left to itself, eventually may decom- 
pose. Just as the energy was able to flow from the 
coordinates of relative motion to some other 
portion of the system, so it may sometimes also 
flow back. If it does this, the energy of relative 
motion assumes its original value, the particles 
are once again free, and they promptly separate. 

This sort of explanation appears quite reason- 
able, but how does one account for the ability of 
the energy to flow in this way? The answer is to 
be found in the type of force existing in a mole- 
cule. It is known that the forces between atoms 
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do not add vectorially and are not central.? The 
potential energy of a group of atoms is a compli- 
cated function of all the coordinates and, strictly 
speaking, cannot be written as a sum of functions 
of individual distances between pairs of atoms. 
Even if one assumes as an approximation that 
the last is possible—so that one can use for the 
individual functions those of the type occurring 
in diatomic molecules, as approximated by 
Morse’ for example—one is still faced with the 
n-body problem. These facts imply that the 
wave equation for a group of atoms is not 
separable for any choice of coordinates. Now if 
the equation were separable in some set of 
coordinates the energy of the system would be 
the sum of a number of energies each due to 
motion involving one coordinate. Each of these 
energies would be constant in time and hence 
there would be no possibility of energy flowing 
from one coordinate to another. Since the 
equation is not separable, however, this is not 
the case. In solving the equation one must use 
approximate methods. One can choose a set of 
coordinates which render the equation nearly 
separable—that is, make it separable if one 
neglects some of the smaller quantities in the 
equation. If one solves this approximate equation 
and then introduces the previously neglected 
terms as perturbations one finds in certain cases 
just the flow of energy described before. This may 
occur, for example, if there are degeneracies in the 
unperturbed problem. 

It is rather a coincidence for a degeneracy to 
occur in the discrete states of a vibrating system, 
that is, for two energy levels belonging to two 
different modes of vibration to coincide, although 
this is known to occur.‘ That discrete levels 
should coincide with continuous is quite common. 
In fact if the system contains two or more 
internal degrees of freedom the energy spectrum 
of each consisting of a group of discrete levels 
surmounted by a region of continuous energy, 
there must be certain discrete energy levels which 
coincide with continuous energies. This can be 


2 London, Probleme der modernen Physik, p. 103, Leipzig 
(1928), Slater, Phys. Rev. 38, 1109 (1931), Dennison, Rev. 
Mod. Phys. 3, 280 (1931). 

3 Morse, Phys. Rev. 34, 57 (1929). 

4 For example in COz, cf. Fermi, Zeits. f. Physik 71, 250 
(1931). . 


seen in Fig. 1, where hypothetical potential 
curves have been drawn for two coordinates so 
that the zero of each is taken as the value of the 
potential energy at an infinite distance from the 
origin. The total energy is the sum of the 
energies for the two coordinates. If now we take a 


Vix, 


Fic. 1. Some energy levels in a system of two degrees ot 
freedom, illustrating degeneracies between discrete ani 
continuous states. 


pair of energy levels Z,,° and E, such that 
their sum is greater than E,°, one can always 
find a continuous level E,,)“? such that 


+ = EnO +E, 


If Z,,° + £, is sufficiently high there may be 
a number of such levels. Thus in Fig. 1, E,.))\! 
and £, represent another such pair. 

In the problem we are considering, since the 
molecule was formed through a collision of two 
free particles, the ‘‘bound”’ or discrete state into 
which it goes (neglecting the possibility of energy 
being radiated, which is slight) must coincide in 
energy with at least one continuous state, namely. 
that in which the particles were before the 
molecule was formed. 

That being the case, there is a certain proba- 
bility for a transition from the discrete state to 2 
continuous one of the same energy—which is 4 
dissociation of the system. This effect in 4 
polyatomic molecule is analogous to the Auger 
effect for an atom or “predissociation’’ for 2 
diatomic molecule, and the method of calculating 
transition probabilities is therefore already avail- 


able.’ The molecule may be said to have 2 


5 Wentzel, Phys. Zeits. 29, 321 (1928). 
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certain lifetime, i.e., the mean time required for 
it to dissociate. If this lifetime is sufficiently long 
so that the molecule can collide with other 
particles before dissociating it may give up its 
excess energy and become a completely stable 
system. 

A case closely related to the one being con- 
sidered at present is the following. There are 
many polyatomic molecules the atoms of which 
are arranged in groups such that the atoms in 
each group are bound together rather tightly but 
the groups are held together quite loosely. By a 
tight binding is ordinarily meant that a plot of 
the potential has a large curvature near its 
minimum. This in general implies that the 
potential valley is deep; whereas a loose binding 
has a shallow minimum and relatively little 
energy is required to cause dissociation. Now if 
such a molecule takes part in a collision as a 
result of which it receives sufficient energy to 
break one of the weak bonds, it may happen that 
this energy is first taken up by one of the strong 
bonds and is present as vibrational energy. After 
some time has elapsed this energy may reach the 
weak bond (in the manner discussed before) and 
cause a decomposition of the molecule. Such an 
effect is of rather frequent occurrence because the 
weak bonds are often in the interior of the 
molecule and thus less exposed to direct effects of 
collisions than some of the strong bonds. This 
process represents a very general type of uni- 
molecular reaction. It is not easy to do much of a 
detailed quantitative nature for this case as yet, 
but perhaps the simple example treated in the 
following section may throw some light on the 
more general problem. 


87? M1 or? 
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where E is the internal energy of the system and 
(4) 


To solve this equation exactly is impossible.® 
One notices however that if one omits the terms 
and the problem reduces 


*For small vibrations one could introduce normal 
Coordinates to get a good approximation. This cannot be 
done in the present case where dissociation is considered. 


II. 
1. General treatment 


We now proceed to treat quantitatively a very 
simple idealized case, our purpose being to 
determine the order of magnitude of the lifetime 
of a molecule having an excess of energy. It is 
hoped that in spite of the simplification the 
essential features of the problem are retained. 

Let us suppose a triatomic molecule with atoms 
of masses m1, m2 and m3, restricted to move along 
a single straight line and let their distances from 
some fixed point on the line be x1, x2 and x; 
(x1 <x2<x3). We make the simplifying assump- 
tion that the potential energy 


U= Vi(x2— + Vo(x3— x2) + (1) 


where, for an equilibrium configuration, | V;| and 
|V2| are large compared to V3. It is to be 
expected that V; and V2 will have forms similar 
to those of stable diatomic-molecule potentials, 
whereas in general V3; will correspond to a 
repulsive force. 

If we write the wave equation for the system, 
introduce new variables 


M1, 

= X3— Xe, (2) 

p= M3X3) /(m+m2+ms), 
and take the wave function as a product ¥(r:, 
r2)x(p), the terms depending on p can be sepa- 
rated off (giving the usual equation for the 


center of gravity of the system) and the equation 
for becomes 


Vi(ri) + Vo(r2) + (3) 


to solving the two ordinary equations 

— (h?/82 /dre+ Vidm® = Em dm”, (5) 
bn /dre?+ Voor (6) 
the solution of (3) being then 


= gm 
(7) 


= + E,®. 
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These quantities will be the more nearly equal to 
the correct ones the larger mz and the smaller V3. 
The operator 


T= Vs, (8) 


which has been omitted can now be introduced as 
a perturbation. 

The quantity V3, as has already been said, is in 
general small. Its effect as a perturbation will be 
to change the energy levels somewhat and it will 
also act to some extent as a coupling between the 
two vibrations characterized by Eqs. (5) and (6). 
However there are few cases in which one knows 
the value of V3. Since furthermore its presence is 
hardly likely to affect the order of magnitude of 
the result we shall not give it much attention, but 
shall center our interest on the first term of 7. 
The effect of this on the energy is zero by a first- 
order calculation if the wave functions are real, 
as is the case here. Its principal réle is as a 
coupling between the vibrations, the physical 
interpretation being that, as the central atom 
moves in both vibrations, it transfers energy 
from one to the other. 

To calculate the rate of transition from a 
discrete state to one in the continuous spectrum 
of the same energy we use the relation® 


T= (40?/h)| SS (9) 


where Wmn° is the wave function for the discrete 
state and y¥.),° for the continuous state 


ds)" denoting a wave function which is a 
solution of Eq. (5) for continuous energy E,.)‘”, 
(i.e., E,.)" >0). The transition rate T' is related 
to the mean lifetime of the molecule + by 


r=1/T. (10) 


2. Numerical calculation 
In order to evaluate (9) we must first get 
solutions of (5) and (6). For V; and V2 we take 
the form suggested by Morse* for diatomic 

molecules 
Vi = Dye — 2D 


| ay) 


V2 = — 2 Doe a2(r2—r02) 


Further we might choose 
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V3 = Dge231+) ( 1 2) 


in accordance with our general ideas about 
interatomic forces. 

In solving (5) and (6) we need our function in a 
somewhat different form from that given by 
Morse and therefore briefly re-derive the solution 
of the equation 


/dr + (82° Dear) 
+2De*"- )Y=0. (13) 


We make the following substitutions: 
(14) 


where we assume that the range of z is 0 to ~ 
(instead of 0 to (42/ah)(2uD)'e"), the error 
being small as Morse pointed out, 


k=(2n/ah)(2uD)}, (15) 
m = 
and (13) becomes 
(16) 


which is the confluent hypergeometric equation.’ 
For large values of z one solution is W,, »(z) 
which has the asymptotic form 


Wr, m(z) {1 (k — 3)? },(17) 


and the other independent solution is W_;, m(—=) 
which obviously does not have the proper 
behavior. We therefore take for our wave func- 
tion »(z). 

We next consider small values of z. The two 
solutions there are M,, »(z) given by 


Mz, m(2) = 
(18) 

and M,, -»(z). If E<0, so that m is real and 

positive, the appropriate solution is Mz, »(2), 

since the other function causes a singularity in ¥ 

at the origin. Now we have the relation’ 

Wa, m(z) —m—k) ]Mz, m(2) 

-m(2), (19) 


7 Whittaker and Watson, Modern Analysis, p. 337, 


Cambridge, 1927. 


r 
t 
a 
| 
fi 
fu 
co 
(1 
pa 


LIFETIMES OF UNSTABLE MOLECULES 


and so our quantization condition is 


=0, 


or 


2tm—k=—n (n=0, 1, 2, ---) 


subject to the restriction m>0. 
We have then for the discrete states 


¥n=2 We, =(—1)"[T(2k —n)/T (2k — 2n) Mz, m(2); 


n nN 
= (2k (—1)i( +A), 
i=0 j 


which is essentially the same as the solution used 
by Morse. 

For E>0, m is imaginary and there is no 
restriction on it. We need the asymptotic form of 
yas or 20. Expanding (19) and keeping 
terms of greatest order of magnitude, one finds 


Ws) ~ 2G cos (sar—srk+y), (22) 
where 


I\(2is)/T(2 —k+is) =Ge'r 
\=arotIn (42(2uD)}/ah). 


23 
and 


(21) 


If one now normalizes the wave function with 
respect to the energy (Appendix, 1, 2) the 
normalizing factor by which it must be multiplied 
is found to be 


= (4/h)(u sinh (24) 
To normalize the discrete wave function, we 


first evaluate a general integral we shall need. 
We find (Appendix, 3) 


n 


f 
j=0 


where y, is a discrete wave-function and y™ an 
arbitrary one.’ For m=k—}—mn and \=0 one 
finds that the normalizing factor for the discrete 
function is 


(26) 


If now one makes use of the easily derivable 
relation 


f Vitn(d V/dr)dr/(En—Ex), (En+Ex) (27) 


and observes that 
dV/dr (28) 


*The notation here is exceptional because the wave 
function is intended to correspond to either a discrete or a 
Continuous state and the superscript is the m defined in 
(15). Everywhere else as far as the appendix, a subscript in 
parentheses indicates a continuous state and is the value of 
5 as defined in (23); without parentheses it indicates a 
discrete state and is the value of the integer in (20). 


G2(n) (2ke—2n 1)/(n!T(2ke —n)). 


(25) 
j'(n—j) (2k 


one finds with the help of (25) that the transition 
rate from a discrete state with approximate wave 
function W»»° to a continuous state with function 
Wn" (if we neglect V3 in (8))° is 


T= (hyde? Fi(s) Fo(m)Gi(p)G2(n), 


(n<m) (29) 
where 


F\(s) = |? sinh 
+ (sinh? rs+ cos? rk;), 
~2|T(kit3+is)|? for s=1, 
F,(m) (2ke—m)(2ke—2m—1), 
Gi(p) = (2k, 


® If one wished to take into account the effect of V; on 
the rate of transition, it could be done by writing (12) in 
terms of 2; and 22 (in analogy with (14)) and then making 
use of the integral (25). This is not done here for the 
reasons previously stated. 
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The quantities s, », p and m are connected by 
the fact that the energy of the two states must be 
the same, which leads to the relation 


(ki — p—2)?) 
= 2)? (31) 
To evaluate (30) we need to know that” 
|? 
MA (32) 


It can be seen from (30) and (32) that if there 
are several continuous states to which the 
system can go, the transition probability is by 
far the greatest to that state which involves the 
least change in the energy of either individual 
vibration (except that if s is very small there is 
a strong dependence on cos 7k;, and predictions 
are difficult, but this is rather unusual). 

Let us now suppose that the molecule was 
formed by the collision of atom 1 with the 
diatomic molecule A consisting of atoms 2 and 3, 


(2k, —2p—1)[T (ki +3) Pm !(2ko— 2m 


each having kinetic energy of the order of 
magnitude of that for room temperature; and let 
us suppose further that as result of the collision 
the vibration of 2 and 3 changed from n=0 to 
some quantum number m, and a vibration of 1 
and 2 with quantum number p was set up. We 
must assume that the original relative kinetic 
energy is such that it matches with the molecular 
levels, i.e., (31) is satisfied. We are not concerned 
with the details of the molecule formation, and if 
we like can assume that a large number of elastic 
collisions take place until the energy match 
occurs. !! 

Then in the decomposition of the molecule the 
process is uniquely defined. We have n=0, s is 
proportional to the initial relative momentum of 
the particles 1 and 2, as in (23), and p is an 
integer satisfying (31). For this case the transi- 
tion rate becomes 


= F\(s)Gi(p) Fo(m) /T(2ke—1). (29a) 


If further we neglect the denominator of (32) 
this can be written 


(29b) 


2m? (2ki — p)T(2ke—1) 


If we set ki p—}3 =¢ we get 


2m —1)[T (ki +3) (2k, —m) 


(29¢) 


In this form we see that I is not as sensitive in 
the dependence of its order of magnitude on k; 
and ke as appeared in the earlier formula. An 
estimate of I may have some reliability (as far as 
order of magnitude) even if there are moderate 
uncertainties in these constants. 

A compound HO, is believed to be formed 
bimolecularly by the reaction of atomic hydrogen 
with molecular oxygen. It is probably not a 
linear molecule and one does not know what the 
constants of its potential functions are. We shall 
assume that the above model can be used 
nevertheless; and for potential functions we use 
those supposed to exist in the lowest electronic 
states of OH and Oz on the basis of their vi- 
brational spectra. It is hoped that the order of 


10 Nielsen, Handbuch der Gammafunktion, p. 23, Leipzig 
(1906). 


magnitude of the result will not be changed very 
much by these approximations. 

We take therefore w2, m2=(16/17, 8, 16) 
1.65-10-*4 g, ky =25.5, ke=68.5, a, =2.0-10° 
cm~!, d2=2.34-108 cm~!. At ordinary tempera- 
tures s is of the order of magnitude of 1. A simple 
calculation leads to the results in Table I. 

We see that for m=1, the calculated lifetime is 
so short that one really cannot speak of the 
molecule having existed, for the entire process 
could not be distinguished from an_ elastic 
collision. For m greater than 4, the lifetime is 
very long compared to the average time between 


1! This matching of energy is not very difficult because of 
the broadening of the levels due to the finite lifetimes of the 
molecular states. Roughly the ratio of favorable to total 
number of collisions is of the order of magnitude of hr/sE 
where AE is the spacing of the vibrational levels in the 
energy region involved. 


Vi 
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collisions at moderate pressures. The mean life- 
time would depend on the collision process which 
would determine the relative frequency of 
occurrence of the various values of m. Because 
collisions (and consequent reactions) would play 
a part in determining the actual lifetime for the 
larger values of m, an experimentally obtained 
value for the mean lifetime would involve a 
certain amount of arbitrariness in its interpre- 
tation. The results calculated here are at any 
rate consistent with available chemical data.'* 


TABLE I. 
m p t I'(sec.~') 
1 21 4 2.4 2-10" 
2 19 6 2-10 
3 17 8 0.5 2-10” 
4 16 9 2.0 4-108 


3. Conclusion 


Aside from the simplification of the conditions 
in the previous calculation, a further difficulty is 


that the unperturbed wave functions do not 
describe the actual situation very well (especially 
for the continuous states) and the perturbation 
may be quite large so that a first order method of 
treatment may be inaccurate. However, one 
would expect the order of magnitude to be 
unchanged in a more accurate treatment. Also it 
must be stressed that the purpose of the calcu- 
lation is mainly to see how the process can come 
about. The particular type of coupling between 
the two unperturbed vibrations does not matter 
very much. The essential point is that in any 
treatment of the problem some sort of coupling 
will be found to exist—and it is to be expected 
that it will in general play a part in bringing 
about the decomposition of the molecule. 

In conclusion the writer wishes to express his 
indebtedness to Professor J. R. Bates, who called 
attention to this type of problem, for discussions 
on the chemical aspects of the subject, and to 
Professor D. M. Dennison for discussions on the 
physical aspects. 


APPENDIX 


1. Normalization of continuous wave functions 


The well-known condition for normalization of a continuous wave function in a one-dimensional problem where the 


variable has a range 0 to is: 


p+A/2 p+A/2 
0 be p—Al2 


where p is a continuous parameter used to order the wave functions. 


Let us write 


Vp =d¢pt+xp, (34) 


where x» is that part of the wave function which goes to zero with sufficient rapidity as x so that fo” xpxpdx 
and fi%¢»x,dx are both convergent. Let us now write J as the sum of two integrals, one for a finite range 0 to a and 
the second for the range a to ~. Then in the first integral and in all the terms of the second which involve x» or xp 


the order of integration can be interchanged and these quantities are seen to vanish. We have therefore 


p+A;2 p+A/2 
J=lim (1/4) | f (35) 
a p—ay2 p—A/2 


where a is finite but can be made as large as one pleases. 
Practically this means that one can use the asymptotic form 
of the wave function in carrying out a normalization. 


2. Normalization with respect to energy 

The normalizing factor for a continuous wave function 
depends on the parameter used. Care must be taken to use 
the same parameter consistently. Thus in the original 


* Bates and Lavin, J. Am. Chem. Soc. 55, 81 (1933); 
Bates and Salley, ibid., p. 110, 


derivation of (9) the energy was the parameter and 
integrations were carried out with respect to it. It is there- 
fore necessary to use the energy in normalizing the function 
in our case. 
If we consider the wave function for a free particle 
er ikz where 
k=P/h=(2uE)'/h, (36) 


P is the momentum of the particle and the range of x is 0 to 
«, the usual normalization with respect to k leads to a 
normalizing factor, 23. If, instead, we use E as parameter 
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the integration is somewhat more complicated and leads 
to a normalizing factor 


N=(1/h)(2u/k)! = (2u/Ph)}. (37) 


For cos 2zkx this is to be multiplied by 2?. Use was made 
of it in (24). 


3. Evaluation of an integral 
To evaluate the integral (25) use is made of Eq. (19). In 


the course of the calculation one obtains two hypergeo- 
metric series for argument 1. These can be summed if 
<1. If 1<A<2, one can make use of the well-known 
relations between hypergeometric functions and by letting 
the argument approach 1, one obtains the same result as 
before. One thus finds that (25) is valid for \<1 and 
1< <2. But since (25) is analytic for 0 <A C2 it must give 
the correct values at \=1 and A=2 also. 
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The lattice energies and lattice constants of the silver and 
thallium halides are calculated assuming ionic crystals with 
a van der Waals potential. The latter term, which is large, 
accounts for the low solubilities (high lattice energies) of 
the salts. The calculations appear to be quantitatively 
satisfactory for the thallium halides, TICI, TIBr and TII, 
and for the three silver halides AgF, AgCl and AgBr. 
Definite evidence is found for assuming the existence of a 
homopolar potential in Ag! of about 10 percent of the total 


lattice energy. One reason for believing this, is the stability 
of the zincblend instead of the rocksalt lattice, which latter 
should be stable were the compound purely ionic. There is 
presumably some homopolar binding in AgBr but it cannot 
be large. The thallium salts probably are entirely ionic. 
The theoretically calculated and experimental (chemical) 
lattice energies are, respectively, in K cal., AgF, 219, 
217.7; AgCl, 203, 205.7; AgBr, 197, 201.8; AglI, 190, 
199.2; TICI, 167, 170.1; TIBr, 164, 165.6; TII, 159, 160.8. 


INTRODUCTION 


HE lattice energy of a uni-univalent ionic 

crystalline salt, namely the energy change 
at O°K of the reaction MX gas, 
can be calculated as the sum of the heats of a 
number of reactions. A series of reactions which 
may be used is shown in Table II. The energy so 
obtained might be termed thé experimental 
lattice energy. If, for the silver and thallium 
halides, these experimental lattice energies are 
compared with the electrostatic energies, it is 
found that the former are almost as large as, or 
even larger than the latter. The electrostatic 
potential is purely attractive, and in addition to 
it there must exist in a crystal another potential 
whose first and second derivatives with respect to 
the lattice constant r (for a uniform increase in 
all crystal dimensions) at the experimental 
equilibrium value of r=7) may be obtained from 
the value of 7) and the compressibility respec- 
tively. 

If this potential of known first and second 
derivatives is assumed to be purely repulsive and 
simple in form (that is the potential V(r)=0 and 
dV/dr=0 for all values of 7, and 9V/dr and 
#V/dr both vary fairly slowly with 1), it is 
impossible to account for the experimental lattice 
energy to within 15 or 20°K cal. for these salts. 
This means that the remaining non-electrostatic 


*A preliminary report was presented at the Congress on 
The Solid State at Leningrad, September, 1932. 


potential is partially attractive, it takes negative 
values at some 7. It is frequently assumed that 
this behavior is to be accounted for by the 
existence of homopolar forces in these crystals. 

There is an additional calculable attractive 
potential besides the electrostatic, the van der 
Waals potential. Its magnitude can be computed 
from equations developed by London! and by 
Margenau.? The calculation used here will follow 
closely a discussion by the author for the alkali 
halides.* If the van der Waals potential is 
included with the electrostatic, the remaining 
potential may be assumed to be purely repulsive 
and the calculated lattice energies agree well for 
all but AgI with the experimental values. The 
method of calculation of the lattice energy 
follows closely the method of Born and Mayer,’ 
whose paper shall be referred to in the succeeding 
paragraphs as B.M. The results of this paper, 
then, show that the homopolar contribution to 
the lattice energies of AgF, AgCl, AgBr, TICI, 
TIBr and TII is practically insignificant, and that 
in AgI is small and can be estimated. 


CALCULATION OF THE VAN DER WAALS POTENTIAL 


The equations used in calculating the potential 
are the following: 


1 F, London, Zeits. f. physik. Chemie B11, 222 (1930). 

2H. Margenau, Phys. Rev. 38, 747 (1931). 

3 J. E. Mayer, J. Chem. Phys. 1, 270 (1933). See April 
number. 

4M. Born and J. E. Mayer, Zeits. f. Physik 75, 1 (1932). 
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Uy = Usat Uga= — (1) 

C= 4 +e _)/2, 

D = Ss’ dy Ss!" (dy _+d_ .)/2, 

2 F 2 F 2 3 P 2 F 2 
3) 


(2) 


3 ete; pi 


Qioo? =— —. (4) 


The complete van der Waals potential in the 
crystal U,, is additively composed of a dipole- 
dipole term Uga varying as r~* and a quadrupole- 
dipole term U,a varying as r~*, with constants 
respectively C and D. C and D can be calculated 
from the corresponding c;;'s and d;;’s, the 
constants for two ions i and j. In order to make 
the calculations the sums S,’, S¢’’, Ss’, Ss’’ must 
be known, which are the sums of /~* and /-® over 
the odd (unlike) and even (like) points of the 
lattice, respectively, if / is the distance between 
the lattice points. The sums depend on the 
lattice type and on the arbitrary choice of 1, 
which here and in the following, is understood to 
be the shortest geometrical constant of the 
lattice, the distance between unlike neighboring 
ions. The sums can be easily computed from a 
table of similar sums given by Lennard-Jones 
and Ingham.® The constants ¢;; and d;; for the 
potential between two ions i and j can be 
calculated from the polarizabilities a and the 
energies ¢ of the ions in the crystal in question. « 
is an energy corresponding to a weighted average 
frequency of the absorption of the ion in the 
crystal.* For the positive ions the polarizability a 
and the energy e are presumed to be practically 
constants of the ions independent of the crystal. 
€ is given the value 90 percent of the ionization 
potential of the gaseous ion (the second ioniza- 
tion potential of the element). For the negative 
ions a is supposed to be inversely proportional to 
e, and ¢ is to be obtained from the experimental 
ultraviolet absorption of the crystal. This as- 
sumption sets Qo”, the 00 matrix element of the 
sum of the g’ for all electrons, a sort of area of the 
ion, as being independent of the crystal. The 


5 J. E. Jones and A. E. Ingham, Proc. Roy. Soc. (London) 
A107, 636 (1925). 
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value of Qo? for the halide ions has been calcu- 
lated from the alkali halides.* », the number of 
outer electrons, is taken as the geometric mean of 
the theoretical (Periodic Table) value, and one 
calculated from the polarizability and frequency. 

For the silver halides the values of ¢ for the 
negative ions have been chosen after an exami- 
nation of the absorption curves for the chloride, 
bromide, and iodide. The value of the polariza- 
bility @ in the crystal has been calculated by the 
use of Qo? and « with Eq. (4). Pauling’s® value 
a=1.72X10-* cm* has been used for Agt and 
the sums of the negative and positive polariza- 
bilities, row 7, Table I, has been compared with 
the experimental value, row 8, calculated from 
the index of refraction m for yellow light. These 
experimental values of }-a (yellow) may be 
expected to be some 5 percent or more higher 
than those which would be calculated were the 
index of refraction extrapolated to X=. An 
estimated value of }-a(A= ©) is given in row 9. 
The agreement indicates only that no great error 
has been made in the choice of constants. 

For the thallium halides the value a=3.50 
X 10-* cm’ for T1* ion is chosen from an exami- 
nation of the index of refraction of some complex 
thallium salts and the corresponding alkali salts. 
The values of ¢ for the negative ions are chosen to 
be in agreement with those of the silver salts. 

The choice of values for all these constants is 
subject to considerable uncertainty. For the 
silver salts the van der Waals potential may be in 
error by as much as 30 percent of its value. For 
the thallium salts the error may be even larger. 
The author has attempted not to err in such a 
way as to make the potential too large. It is felt 
that the actual potential is more likely to be 
greater than to be smaller than the values 
calculated here. Owing to the use of an empirical 
repulsive potential, determined from the at- 
tractive potential and the constants of the 
crystal, the magnitude of the error introduced 
into the computation of the lattice energy by an 
error in the assumed values of C, and especially 
of D, is less than that introduced into the van der 
Waals potential. That is, if the value of Uw 
= —C/r were assumed to be 3 k cal. greater 
than the values here given, say for AgF, the total 


* 1. Pauling, Proc. Roy. Soc. (London) A114, 181 (1927). 
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calculated lattice energy would be only 1.5 k cal 
greater. 


THE REPULSIVE POTENTIAL AND THE 
LATTICE ENERGY 


The compressibilities, 8, of AgCl, AgBr, 
AgI, TICI and TIBr are given in International 
Critical Tables. The thermal expansion for 
AgCl and AgBr is also given there. The term 
(37/ NB)(dv/T), in which the thermal expansion 
is needed, is small, arid its value can be estimated 
for the other salts without introducing any 
appreciable error. From these data and the 
lattice constant, 7, which is taken from the 
Strukturbericht’ the repulsive potential is com- 
puted. The method deviates from that of B.M. 
only in so far that for the silver salts it is not 
assumed that the exponents for the repulsive 
potential between two Agt ions and two halide 
ions, or an Agt ion and a halide ion, are the 
same. 

If the repulsive potential energy of the crystal 
as a function of 7 is written B(r) it has been 
shown that, at r=”, the lattice constant at 
temperature 7, and P=0, 


dB(r) ae? 6C 8D 3T saV 
dr Yo ro° 
42C 72D 9V 


dr* Yo NB 


(6) 


These equations differ from those in B.M. only in 
the inclusion of the D-term and the omission of a 
factor multiplying 9V//N@ in Eq. (6). This factor 
always has a value very close to 1. 

It is now assumed that the replusive potential 
between two ions i and j, b;;(r;;) is given by the 
equation, 


bij = je"! 4, (7) 
with 
= Pi +p;, 
2pu- = p_ = 0.345 X 10-8 cm, 
(8) 
2pagt = p+, 
pu-+pagt =p. 


a’ P. Ewald and C. Hermann, Strukturbericht, 1913- 
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Eq. (7) is that used in B.M. and Eq. (8) is 
introduced to allow for the fact that not all ion 
pairs show the same value p;;, which, however, 
was the case for the alkali halides. There is some 
reason to believe that Eq. (8) is justified. It is 
sufficient here to point out that the same 
repulsion is assumed between two negative ions 
that was found in the alkali halides. The re- 
pulsion between the positive ion pairs is negligible 
for all silver salts on this, or any other reasonable 
assumption, and the value of p, is in reality never 
used. The value of p, the repulsive exponent 
between an Agt and a halide ion, is calculated 
from the compressibility, and the assumptions 
used are relatively unimportant in influencing the 
value obtained. That is, actually Eq. (8) is 
entirely unimportant, the repulsion between Agt 
pairs is probably negligible and is neglected, that 
between halide pairs is assumed to be the same 
function of 7 as in the alkali halides, and that 
between Ag+ and halide ions is calculated from 
the constants of the crystal. 

Defining as M the number of unlike neighbors 
of any lattice point, as M’ the number of like 
“‘next-to-neighbors,” as a the ratio of the 
distance between next-to-neighbors to that be- 
tween neighbors, and with ¢,;=¢,4/c,_, Ce 
=c__/c,_, the relative effect of the charge on the 
repulsion, so that c,;=1.20 and ce=0.81 and 
c, _=0.929 is written co, one obtains for the total 
repulsive potential in the crystal, 


B(r) = beol + 
(9) 


M,= M. bcobs.b_e~!?, 


2Mb_ 
M "Cob_ 


The constant a has the value 1.155, 1.414 and 
1.633 for the CsCl, NaCl and Zincblend lattice, 
respectively. Co, ¢; and C2, the dependence of the 
repulsion on the charge are computed from a 


M2 


theoretical equation of Pauling’s* assuming 18 


8 L. Pauling, Zeits. f. Krist. 67, 377 (1928). 
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electron shell positive ions and 8 electron shell 
negative ions. The values of the b_’s are calcu- 
lated from the alkali halides, and those of the 
b,’s are later calculated in this paper. The 
constant 0 is arbitrarily chosen to be 10-” erg. 

M,* represents the total repulsion due to pairs 
of Ag* ions divided by that due to Agt halide 
pairs, and M,.~ the total repulsion of halide pairs 
divided by that of the Agt halide pairs. M.2* at 
r=ro is entirely negligible for all silver salts and 
will be omitted from the further calculation. 

For the thallium salts p is found to have very 
nearly the same value as in the alkali halides and 
the somewhat simpler equations of B.M. have 
been used throughout. The equations are equiva- 
lent to Eq. (9) setting p=p_=p, and adding 
M ot+ 2 = M 2. 

Using Eq. (9) one obtains for the derivatives, 


Yo 
=>— (10) 


T=-—71% 


or p 
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TABLE I. Calculation of the theoretical lattice energies. 


with 

y=ap/p_, 
when MM, and M-.~ are understood to be the 
values of these functions at r=”. This leads to 
the equation for 


ro/o= (0/7) Mar (13) 


Table I includes many of the numerical steps 
of the calculations and the value of p, row 33. 
An average value, p=0.260 10-8 cm and 0.361 
x<10-% cm for the silver and thallium salts, 
respectively, is obtained. The further calculations 
are made with these average values. 

B(r.) the repulsive potential at r=7) may now 
be calculated from Eqs. (9) and (10), 


B(ro) =(1—k)r(0/10), (14) 


1 Salt AgF AgCl 
2 Crystal ty NaCl NaCl 
3 Qoo? (neg. ion) X 10" cm? 1.44 3.40 
4 e« (neg. ion) X10” ergs/molecule 19.0 15.6 
5 @ (neg. ion) X 10% cm? 1.15 3.45 
6 @ (pos. ion) X 10% cm? 1.72 1.72 
7 La (calc.) cm* 2.87 5.17 
8 Za (yellow) 10% (exp.) cm? — 5.41 
9 Za (A= ©) X10* (exp.) cm® 5.14 
10 (pos. ion) X 10" ergs/molecule 30.0 
11 4X10 ergs 67 — 
“ 19 133 
13 c, “ 34.5 89 
14 d, 10% ergs cm® 91 — 
23 268 
16 d, * 44 150 


21 CX 10 ergs Xcm® 
22 DX10" ergs X cm’ 


23 ro X 108 - em 2.46 2.772 
24 8X10° (compressibility) bars™ — 2.40 
25 V-\(0V/dT)p X 108 — 98.9 
26 9V/NB X10” ergs/molecule — 159 


27 (37 /NB)(0V/dT) X10" ergs/molecule (1.50) 1.57 
28 ae?/ro X10” ergs/ molecule 16.20 14.39 


29 C/ro&X102 / 1.39 1.70 
30 D/r&x10% “ / “ 0.24 0.32 
31 “ / 24.9 25.6 
32 * — 283. 
33 a/t 11.3 
34 pX108 cm 0.247 
35 ro/p calculated 9.45 10.65 
36 B(r) X10” ergs/molecule 2.63 2.38 
37 0.10 0.08 
38 U X10" ergs/molecule 15.10 13.95 
39 U k cal. 219 203 


AgBr Agl TICI TIBr TI 
NaCl = Zincblend CsCl CsCl CsCl 
4.08 $.27 3.40 4.08 5.27 
13.9 11.0 15.6 13.9 11.0 
4.47 7.29 3.45 4.47 7.29 
1.72 1.72 3.50 3.50 3.50 
6.19 9.01 6.95 7.97 10.79 
— — 268 
208 437 133 208 437 
109 151 179 222 307 
475 1228 268 475 1228 


2.884 2.812 3.33 3.44 3.62 
2.75 4.11 4.9 5.3 — 
156 150 104 106 —— 
1.63 (1.70) (1.60) (1.70) (1.80) 
13.82 13.29 12.04 11.66 11.09 
1.60 1.74 1.66 1.67 1.75 
0.29 0.38 0.24 0.25 0.34 
24.6 25.1 22.3 22.0 22.5 
276 277. 215. 217. — 
10.8 11.0 9.64 9.90 — 
0.271 0.260 0.356 0.366 — 
11.09 10.78 9.21 9.51 10.03 
2.20 2.28 2.36 2.23 2.14 
0.06 0.05 0.06 0.05 0.04 
13.53 13.08 11.52 11.30 11.00 
197 190 167 164 159 
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TaBLe II. Calculation of the experimental lattice energies. AH in kilocalories of the reactions.* 


Reaction AgF AgCl AgBr Agl TICI TIBr TU 
MXer, 0°K—M Xer, 291°K 2.5 2.8 2.9 3.0 2.9 3.1 3.2 
MXer.~Ma.s.+3X25.s. 291°K 48.7 30.6 23.7 14.9 48.7 41.1 30.1 
3X2 5.8.7 Xz. 291°K 31.8 28.8 26.7 25.5 28.8 26.7 25.5 
Xg. 291°K—Xg. 0°K —1.1 -1.1 —1.1 —1.1 —1.1 —1.1 —1.1 
Xp. te7—>Xg.~ O°K —95.3 —86.5 —81.5 —74.2 —86.5 —81.5 —74.2 
Mer, 291°K—Mer, 0°K —1.7 —1.7 —1.7 —1.7 —1.7 —1.7 —1.7 
Mer. Mg. 0°K 59.0 59.0 59.0 59.0 39.0 39.0 39.0 
Mg.>Mg.*+e~ 0°K 173.8 173.8 173.8 173.8 140.0 140.0 — 140.0 
Sum = M Xer.—>Mg.+ + O°K 217.7 205.7 201.8 199.2 170.1 165.6 160.8 
Theoretical lattice energy 219 203 197 190 167. 164. 159. 
Difference —1.3 2.7 4.8 9.2 2.9 1.6 Be 


* Subscripts, cr. indicate crystal; s.s., standard state; and g, gas. 


The complete lattice energy U is given by 


ae C D 
-—-—- —+B(ro)+3hv, (15) 
the small value of hy)/2 can be easily estimated. 
The values of the lattice energy in ergs per 
molecule and in k cal. are tabulated in the last 
two rows of Table I. 

In Table II are given the energies of the 
reactions used to calculate what is called the 
experimental lattice energies. The heats of for- 
mation of the crystals are taken from Jnter- 
national Critical Tables, as are also the heats of 
formation of the halogen gaseous atoms. The 
small energy changes between 291°K and 0°K are 
all estimated. The electron affinity of the halogen 
atoms is taken from a paper by Mayer and 
Helmholz® where they were calculated from the 
lattice energies of the alkali halides. The heat of 
sublimation of the metals has been estimated by 
the use of the Hildebrand” modification of 
Trouton’s rule. The sum of the AH’s of all of the 
first eight reactions give the AH, which at 0°K is 
the AE, of the dissociation of the crystal into 
gaseous ions, the lattice energy. The theoretical 
lattice energy and the difference between the two 
values is also tabulated. 
The uncertainties in the lattice constants and 
in the compressibilities may account for 1 k cal. 
uncertainty in the theoretical lattice energies. 
The van der Waals potentials of about 20 k cal. 
may be in error by as much as 6 k cal., which 
would contribute about 3 k cal. error to the 
lattice energy calculations. In all the values of the 
i Mayer and L. Helmholz, Zeits. f. Physik 75, 16 

J. H. Hildebrand, J. Am. Chem. Soc. 37, 970 (1915). 


theoretical lattice energy may be in error by 
about 4 k cal. as an upper limit, even if the 
assumptions under which they are calculated 
hold exactly. The experimental lattice energies 
may also be in error by as much as 3 or 4 k cal., 
but the difference between the values for different 
halides of the same positive ion are not likely to 
be in error by more than 2 k cal. The trend of 6k 
cal. in the difference between the theoretical and 
experimental energies from fluoride to bromide of 
the silver salts is probably real, and represents a 
deviation from a purely ionic character. The 
trend could be diminished by assuming a greater 
polarizability for the Agt ion, or in general by 
changes in the assumptions tending to increase 
the van der Waals potential; but it seems 
probable that this trend shows an increasing 
homopolar contribution to the lattice energy 
which is zero for the fluoride and about 2 or 3 
percent of the total energy in the bromide. 

This view is strengthened by the marked 
difference between the theoretical and experi- 
mental energies in the silver iodide. It seems that 
there is no doubt that this difference cannot be 
explained by errors in the choice of the van der 
Waals potential. Further evidence will be de- 
duced later in this paper to show that the 
homopolar bond contributes approximately 10 
percent of the total experimental lattice energy. 

The thallium halides on the other hand, are 
seen to be all, within the accuracy of these 
calculations, purely ionic. 


Tue LATTICE CONSTANT 


The Eggs. (7) and (8) assume that the repulsive 
potential between any two ions 7, j, can be 
calculated from constants of the ions, );, b;, pi, 
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TABLE III. Values of theoretical lattice constants. 


TIBr TI 


AgCl AgBr Agl TICI 
Value of b_ 26.80 77.4 105.3 189.5 774 105.3 189.5 
Value of b, 187.8 187.8 187.8 187.8 33.78 33.78 33.78 
Calculated ro X 108 cm 2.45 2.78 2.89 2.952 3.30 3.42 3.69 
Experimental ro X 108 cm 2.46 2.772 2.884 2.812 3.33 3.44 3.62 
Difference, 67 X 108 cm —0.01 0.01 0.01 0.14 —0.03 —0.02 +0.07 


p;, and c;; the dependence on charge and outer 
electron configuration. Actually the method used 
here and in B.M. does not do this. The constant p 
is indeed calculated from several crystals and an 
average value chosen for the repulsive potential, 
but the product 0,; is implicitly determined for 
each crystal by the use of the lattice constant. 
‘The ratio b;/b; does enter the equations used, but 
in such an unimportant way (in the computation 
of the small 1,+ and M,~) that in the calculation 
of the lattice energies of the alkali halides the 
Goldschmidt radii were used to estimate the 
ratio, and no attempt was made, or needed, to 
justify this, except in a most qualitative way. 

Huggins and Mayer!! have shown that 9 
values of b for the 9 ions of the alkalis and halides 
can be so chosen that one can calculate the 20 
lattice constants of the 20 alkali halides to 
within about 0.01A, probably within the experi- 
mental error of their determination. 

That is, Eq. (7) is justified experimentally for 
the purely ionic alkali halides. With the values of 
b_ for the halides found from the alkali halides, 
the values of 6, for Ag+ and Tl*+ may be calcu- 
lated. One will obtain 4 values from the 4 halides 
for Ag+, and 3 values from the 3 halides for TI*. 
The values of b for TI* are found to have a slight 
trend, which is probably of significance but will 
be ignored, and the average of all used to obtain 5 
for Tl*+. The values of the 6 for Ag* decrease 
noticeably with increasing atomic number of the 
negative ion. Assuming that the salts AgF and 
AgCl are essentially ionic, the true b, for Agt 
may be taken as an average of that determined 
from these two salts. (The lattice constant of AgF 
is sufficiently uncertain to warrant the use of 
AgCl although it may be somewhat homopolar.) 

With these average values of b for Ag* and TI* 
the lattice constants of the salts may be calcu- 


11M. L. Huggins and J. E. Mayer, Phys. Rev. 43, 373A 
(1933). Abstract of Report at Atlantic City Meeting, 
December, 1932. 


lated. The results are tabulated in Table III. 
For the thallium salts although the calculated 
lattice constants do not agree well with the 
experimental values, the differences are probably 
not much greater than the experimental un- 
certainty in the lattice constants, which are not 
well known. A trend is apparent in the differences. 
For the silver salts the trend is decided and large. 
The calculated theoretical lattice constant for 
Ag] is considerably higher than the experimental. 

This small experimental lattice constant is in 
good agreement with the observations on the 
lattice energy for AgI. A homopolar bonding 
potential is probably pulling the ions closer 
together and simultaneously lowering the energy. 
We may regard the calculated lattice constant as 
the lattice constant of a purely ionic crystal. The 
lower experimental lattice constant is due to an 
actual existent homopolar binding in addition to 
the ionic. 

One may now compute more exactly the 
homopolar part of the potential in AgI. The 
repulsive potential in Table I has been calculated 
from the empirical lattice constant of the crystal. 
Using instead the value of 6 for Ag+ found from 
the fluoride and chloride one may calculate the 
repulsion in a truly ionic crystal. This truly ionic 
crystal would have the lattice energy 184 k cal. at 
its equilibrium lattice constant 2.95A, and 178 k 
cal. at the experimental lattice constant 2.814, 
corresponding to the highly compressed crystal. 
Subtracting this from the experimental lattice 
energy one arrives at a value 21 k cal., or 10 
percent, which is due to the existence of a 
homopolar potential. 

The actual calculations of the 5,’s and the 
lattice constants are not reproduced here in 
detail. By using Eq. (9) and the value of B(r) at 
r =r given in Table I, b, can be calculated if b_ is 
known. Having an average value of 5, the value 
of B(r) at r=r can be recalculated. In order to 
calculate the theoretical lattice constant with the 
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new theoretical value of the repulsion B(r) one 
seeks a value of r that will satisfy Eq. (5). By 
using a method of variations the following 
approximation is obtained. If bo is the average 5 
which we have chosen as the correct b of the 
positive ion, and 0}, the experimental value 
previously determined from the lattice constant 
of the crystal, one finds for 6r, the difference 
between the theoretical lattice constant and the 
experimental: 


dr = p(w—1)/(9V/NBo) (16) 


with 
o= bo/de. 


THE STABILITY OF DIFFERENT LATTICE TYPES 


In the previous discussion it has been assumed 
that insofar as the compounds studied are purely 
ionic and not in such a crystal that the ions are 
not deformed (which means that one is restricted 
to regular crystals), the complete potential be- 
tween any two ions is known as a function of the 
distance between them. In principle one can then 
calculate the lattice energy of any regular crystal 
lattice for these compounds. The actual work of 
solving for the equilibrium lattice constant 7 out 
of Eq. (5) would be somewhat tedious. 

It is, however, comparatively simple to obtain 
the change of energy with change of lattice type 
at constant 7, and the change of energy with r of 
the new hypothetical lattice, as it adiabatically 
expands or contracts to its equilibrium 7 can be 
easily estimated. The method is that used on the 
alkali halides in the paper B.M. For each term, 
except the hv)/2, of the expression for the lattice 
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energy, Eq. (15), one may calculate the change in 
the value of that term due only to a change in the 
lattice type, but with 7) held constant. The sum 
of these changes will give the complete variation 
(6@),.,, in the lattice potential with change of 
crystal type at constant 7. The small term, 
invariable negative, corresponding to the de- 
crease in lattice potential, }(0°/dr*)(6r)?, as the 
new lattice type changes to its equilibrium 
density, is then added. Starting always from the 
experimental type a net positive value for the 
total 6 indicates that the calculation shows the 
hypothetical lattice type to be instable. A net 
negative 6 indicates an error of some sort in the 
calculation, provided the hypothetical lattice 
type is not actually experimentally obtained. In 
Table IV the results of such a calculation are 
shown. 

It is to be seen that AgF, AgCl and AgBr are 
calculated to be instable as Zincblend or CsCl 
lattices, and should, as they do, crystallize in 
rocksalt lattices. AgI, on the other hand is 
calculated to be stable in the rocksalt lattice, 
which is in contradiction to the experimental 
evidence. The discrepancy of 0.43 X 10-" erg or 6 
k cal. is far too large to be accounted for by 
anything but a failure of the fundamental 
assumption that it is ionic. The salts TICI and 
TIBr are found to be stable as rocksalt lattices, 
although they actually crystallize in the CsCl 
lattice. Here, however, the discrepancy of 0.02 
and 0.01X10-™ erg is so small that it is easily 
accounted for by a relatively small error in the 
van der Waals potential. TII is found to be stable 
in the CsCl type lattice. 


TABLE IV. Calculation of potential energy change with change of lattice type. 
Energies in 10~ erg/molecule. 


Experimentally Hypothetical 


Change at constant 7 of 


stable lattice lattice 


—ae?/ro 


—C/r& | 


Rocksalt Zincblend 


Zincblend Rocksalt 
Rocksalt CsC 


“a 


CsCl Rocksalt 


+1.02 
+0.90 
+0.87 
—0.89 
—0.14 
—0.12 
—0.12 
+0.10 
+0.10 
+0.09 


+0.55 
+0.67 
+0.67 
—1.23 
—0.67 
—0.82 
—0.81 
+0.56 
+0.55 
+0.57 


+0.09 
+0.13 
+0.11 
—0.26 
—0.12 
—0.16 
—0.15 
+0.09 
+0.08 
+0.11 
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AgF | -0.95 | —013 | +0.58 
AgCl —0.96 | -013 | +0.61 
in AgBr —0.93 —0.13 +0.49 
Agl +2.08 —0.13 —0.43 
A AgF +1.16 —0.06 +0.17 
is AgCl +1.72 —0.06 +0.56 : 
1c AgBr +1.99 —0.06 +0.85 
TICI —0.73 —0.04 —0.02 
0 TIBr —0.70 —0.04 —0.01 
TH —0.70 —0.04 +0.03 


DISCUSSION OF THE RESULTS 


The salts here considered all show two distinct 
variations from the characteristics of the defi- 
nitely ionic alkali halides. Their lattice energies 
are markedly high, which is qualitatively evi- 
denced by their low solubilities, and their lattice 
constants do not increase as rapidly with in- 
creasing halide atomic number as do the alkalis. 
The difference in lattice constant of AgF and 
AgBr is 0.42A as compared with 0.67A for the 
difference between NaF and NaBr. The difference 
ro(TIC1) —7o(TII) is 0.29A. Compared with 0.39A 
for CsCl and CsI. 

A high attractive van der Waals potential, 
increasing with the size of the negative ion, and in 
the case of the silver salts a low repulsive 
exponent p deducable from their low com- 
pressibilities, qualitatively explain both of these 
anomalous characteristics. In the case of the 
three salts AgF, AgCl and AgBr the explanation 
is satisfactory quantitatively. There is some 
evidence that AgBr is not entirely ionic, as one 
would expect if AgI is supposed to have an 
appreciable homopolar potential, but it would be 
possible to account for the deviations from the 
calculations by assuming a somewhat greater van 
der Waals potential. AgI, however, cannot be 
accounted for quantitatively as ionic. Not only 
are the deviations from the lattice energy and 
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constant too great to be patched by a small 
alteration of the assumed polarizabilities and 
frequencies, but more important, it crystallizes in 
the wrong type lattice. Were it ionic, and 
especially if the van der Waals potential were 
larger, it would be instable in the zincblend 
lattice, (and also in the similar Wurtzite lattice). 
The homopolar potential apparently favors 
markedly the low coordination number type of 
lattice, which is probably to be expected. 

The quantitative calculation of the TICI, TIBr 
and TII lattice energies is satisfactory. The 
anomalous lattice constants are not accounted 
for successfully, although the experimental values 
are sufficiently uncertain to allow some latitude. 
The calculated lattice constants would agree 
better with the experimental, as would also the 
lattice energies, were the van der Waals potential 
assumed to be greater, particularly for the iodide. 
Here the study of the stability of different types 
corroborates the idea that the van der Waals 
potential is higher than that used. A higher van 
der Waals potential is needed to account for the 
experimental stability of the CsCl type lattice, 
whereas an existent homopolar potential would 
be expected to favor the lower coordination 
number of the rocksalt lattice. The author con- 
siders this a strong reason for believing that prac- 
tically no homopolar binding is present in these 
salts. 
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Some Conclusions from Rules on the Concentration of Valency Electrons in Binary 
Intermetallic Alloys 


HARALD PERLITz, University of Estonian Republic 
(Received January 12, 1933) 


1. INTRODUCTION 


N 1926 W. Hume-Rothery suggested that the 

similarity of B-phases in Cu-Zn, Ag-—Zn, 
Au-Zn, Ag-—Cd, Au-Cd, Cu-Al, Ag—Al and 
Cu-Sn alloys is due to the fact that these phases 
are space lattices of atoms and electrons for 
which the ratio of electrons to atoms, the so- 
called concentration of valency electrons, is the 
same, viz., 3 : 2, according to the chemical for- 
mulae CuZn, AgZn, AuZn, AgCd, AuCd, Cu;Al, 
Ag;Al and Cu;Sn which are usually ascribed to 
these 6-phases.! At the same time, from a survey 
of displacements of homogeneity ranges of struc- 
turally similar phases in certain binary alloys of 
copper, silver and gold Arne Westgren and Gésta 
Phragmén concluded too that the atomic group- 
ing of intermetallic phases is to a considerable 
degree regulated by their concentration of 
valency electrons.’ 


2. STATEMENT OF THE RULES 


To-day 35 binary intermetallic alloys are 
known in which up to three intermediate phases 
are met with which are structurally similar to 
B-, y- and e-phases of the Cu-Zn-system.’ As a 
rule these phases occur at such compositions for 
which the concentrations of valency electrons 


‘J. Inst. Met. 35, 295 (1926), London. 

* Ark. mat., astron. och fysik B19, 5 (1926). 

* Arne Westgren und Gésta Phragmén, Zeits. f. Metallkde. 
15, September (1926); Metallwirtsch. 7, 701-702 (1928); 
Zeits. f. anorg. allg. Chemie 175, 84 (1928); Zeits. f. 
Metalikde. 18, 279 (1929); Trans. Faraday Soc. 25, 380- 
382 (1929). A. Westgren und W. Eckman, Ark. kem., 
mineral. och geol. B10, No. 11, 1-4 (1930). W. Eckman, 
Zeits. {. physik. Chemie B12, 58-59, 77 (1931). A. West- 
gren, Zeits. f. Metallkde., No. 11 (1930); Metallwirtsch. 9, 
922-923 (1930); Trans. Amer. Inst. min. metallurg. Engr. 
Institute of Metals Div. 1931 p. 35-37; J. Frank. Inst. 212, 
592-595 (1932); Zeits. f. angew. Chem. 45, 11, 13, 15-18 
(1932); Harald Perlitz, Acta et Commen. Univ. Tartuensis, 
24, tract No. 3 (1932); Metallwirtsch. 12 (1933). 


attain certain values peculiar to each type of 
structure. The characteristic values for concen- 
trations of valency electrons are: 

(1) 3:2 for phases with the body-centered 
cubic lattice of 8-brass or with the cubic lattice of 
8-manganese with 20 atoms in the unit cube, 
phases which will be referred to as 8-phases. 

(2) 21 : 13 for phases with the cubic lattice of 
y-brass with 52 atoms in the unit cube, phases 
which will be referred to as y-phases. 

(3) 7:4 for phases with the close-packed 
hexagonal lattice of e-brass, phases which will be 
referred to as e-phases. 


3. FORMULAE FOR AND FROM THE RULES 


Now, if the number of valency electrons of 
the first component of the binary alloy is denoted 
by m, and that of the second one by 2; if in the 
phases with B-, and y-, and e-structures the 
number of atoms of the first component is 
denoted by x, and that of the second one by y, 
then the concentration of valency electrons, , is 
defined by (mx+ny) : (x+y)=k. Solving this 
equation with respect to the ratio x/y, we have 
x :y=(k—n) :(m—k). In this expression the 


‘left-hand side is positive, and therefore the right- 


hand side should be positive too. But that 
necessitates that m, n, and k should satisfy the 
inequality mSk=n. 


4. CONCLUSIONS FROM THE INEQUALITY 


The obtained inequality leads to following 
conclusions: 

(1) B-, y-, and e-structures are not to be 
expected in binary intermetallic alloys if both 
components have the same number of valency 
electrons. 

(2) AB-, y- or e-structure is to be expected in 
binary intermetallic alloys if the one component 
has less and the other component more than 
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3:2, 21:13 or 7:4 valency electrons, respec- 
tively. 

(3) Since the characteristic values of concen- 
trations of valency electrons are between two 
consecutive integers, it follows that in alloys in 
which one of the three structures referred to is 
met with, both other structures should be ex- 
pected too. This statement accounts for the 
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frequent coexistence of B-, y-, and e-structures. 

(4) Since the lower and upper limit for the 
values of concentration of valency electrons are 
one and two, it follows that a 8-, y- or e-structure 
is to be expected in binary intermetallic alloys 
if the one component has not more than one 
valency electron, and the other component has 
not less than two valency electrons. 
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two mercaptans are measured and compared with values in 


the literature. The differences in moment of corresponding 


oxygen and sulfur compounds are attributed to inherent 


(Received February 28, 1933) 


The dipole moments of thianthrene, two thioethers and 


differences in the H-O, C-O, H-S and C-S moments and 
to the greater polarizability of sulfur as compared to 
oxygen. In contrast to the oxygen ethers, the thioethers 
show as much intermolecular action as the mercaptans. 


HIS investigation was undertaken some 

time ago with the object of comparing the 
mercaptans and sulfides, for certain of which 
approximate moments had been calculated pre- 
viously, with the corresponding oxygen com- 
pounds, which had shown unusual dielectric 
behavior.2, Hunter and Partington, who de- 
termined the dipole moments of ethyl mercaptan 
and ethyl sulfide some little time ago,* have just 
published values for thiophenol, two alkyl 
mercaptans, and three more sulfides,‘ and have 
given an interesting discussion of their results. 
Although their values for the sulfides are in 
excellent agreement with those in the present 
paper, as will shortly be seen, there is a marked 
discrepancy in the values for the mercaptans, 
which are difficult to obtain in a very pure state. 
On this account an unusually detailed description 
of the criteria of purity of the two mercaptans 
here measured seems desirable. 


PURIFICATION OF MATERIALS 
Benzene 


Mallinckrodt’s (Thiophene Free) analytical 
reagent grade benzene was treated as in earlier 
work.® 


'C. P. Smyth, J. Am. Chem. Soc. 46, 2151 (1924). 
*C. P. Smyth and W. N. Stoops, J. Am. Chem. Soc. 51, 
3312, 3330 (1929). 


*E. C. E. Hunter and J. R. Partington, J. Chem. Soc. 
1931, 2062. 


*E. C. E. Hunter and J. R. Partington, J. Chem. Soc. 
1932, 2812, 2819, 


°C. P. Smyth and W. S. Walls, J. Am. Chem. Soc. 54, 
1854 (1932), 
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Thianthrene 


Material very kindly loaned by Professor 
Gregg Dougherty (m.p. 156.5-157.5°) was frac- 
tionally crystallized nine times from pure dry 
benzene and dried in an evacuated desiccator. 
This material was odorless and concentrated 
solutions in benzene showed no coioration; m.p., 
157.5-157.8°. Bergmann and Tschudnowsky re- 
port® their material to melt at 154°. 


Ethyl sulfide 


Material from the Eastman Kodak Company 
was extracted with sodium hydroxide solution, 
washed several times with water, shaken with 
mercury (no decoloration), dried for several days 
with anhydrous sodium sulfate, decanted and 
twice fractionally distilled; b.p., 91.8° (756 mm) ; 
np” 1.44265; np* 1.4400. 


n-Amyl sulfide 


Material from the Eastman Kodak Company 
was treated in the same way as the ethyl sulfide 
and twice fractionally distilled under reduced 
pressure; b.p., 103—103.5° (13 mm); mp*° 1.45547; 
np” 1.45377. 


n-Butyl mercaptan 

Material from the Eastman Kodak Company 
was dried for several days with anhydrous sodium 
sulfate and fractionally distilled in an apparatus 
filled with nitrogen. The boiling point was con- 
stant at 97.5° (760 mm) during the collection of 
the first three fractions, after which it rose 0.1°. 
A 240 cc sample was collected in seven fractions 


6E. Bergmann and M. Tschudnowsky, Ber. 65, 457 
(1932). 
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in the quantities given below with the following 
refractive indices: (1) 50 cc, 1.4381; (2) 40 cc, 
Np” 1.4381, np” 1.4356; (3) 35 cc, mp” 1.4381; (4) 
50 cc, np?” 1.4387; (5) 25 cc, mp” 1.4402; (6) 20 cc, 
np” 1.4416. The second fraction, which was used 
for preliminary dielectric constant and density 
measurements, had an index considerably lower 
than the value 1.4411 found by Hunter and 
Partington and a density, 0.8377, obtained by 
extrapolation to 20°, also lower than the value 
0.8397 found by Hunter and Partington and the 
value 0.8412 obtained from the specific volume 
determined by Bingham and Fornwalt.’? The 
density of this second fraction was, however, 
slightly higher than the value 0.836 given in 
International Critical Tables and its boiling point 
97.5° was between the value 98° given by 
International Critical Tables and the value 97.2° 
found by Hunter and Partington. 

In order to remove possible butyl alcohol, 
fractions 1 and 4 combined were washed five 
times with distilled water, dried with anhydrous 
sodium sulfate, further dried with phosphorus 
pentoxide for a few hours, decanted, and frac- 
tionally distilled as before; b.p., 97.9° (755 mm); 
Np” 1.44275; np™ 1.44025; d,?° 0.8415. These 
constants are in excellent agreement with those 
recently published by Ellis and Reid.® 

As this sample appeared to be purer than the 
first material measured, it was used to obtain the 
data given in this paper. The dipole moment 
obtained for it was 0.04 10~'* lower than that 
found for the first sample. 


n-Amyl mercaptan 

Material from the Eastman Kodak Company 
was dried for several days with anhydrous 
sodium sulfate, decanted, and fractionally dis- 
tilled three times; b.p., 123.7-124.0° (756 mm); 
np” 1.4459; nmp* 1.44386. The density, 0.8418, 
obtained by extrapolation to 20°, is in excellent 
agreement with the value 0.8415 of Bingham and 
Fornwalt, but much lower than the value 0.857 in 
International Critical Tables, which gives a some- 
what lower refractive index, 1.4437, and a higher 
boiling point 126°. Our values are in agreement 


7E. C. Bingham and H. J. Fornwalt, J. Rheology 1, 372 
(1930). 

8 L. M. Ellis, Jr., and E. E. Reid, J. Am. Chem. Soc. 54, 
1674 (1932). 
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with the values d,;” 0.8375 and mp” 1.4440 
obtained by Ellis and Reid. 


EXPERIMENTAL RESULTS 


The densities d and dielectric constants € of the 
polar substances of mole fraction cz in solution in 
benzene were measured with the apparatus and 
methods previously employed,’ a wave-length of 
1000 meters being used in the dielectric constant 
determinations. The results are shown in Table I, 
in which the values opposite ce=0 for the ben- 
zene-thianthrene mixtures are for the pure 
benzene used in all the measurements. P, 
=(e—1)M/(e+2)d has been extrapolated to 
c2=0 in the usual manner to obtain P., from 
which, together with the molar refraction 
MRp, the moment is calculated as h=0.0127 
X10-8[(P,.— MRp)T]}!. The values are sum- 
marized in Table II, and the moments of the 
mercaptans and sulfides are shown in Table III 
marked ‘“‘W. and S.”’ together with the values of 
Hunter and Partington in the column marked 
“H. and P.” 


DISCUSSION OF RESULTS 


The moment, 1.41X10-'8, found for thi- 
anthrene is decidedly lower than the value, 
1.68 <10-"8, just published by Bergmann and 
Tschudnowsky’ but surprisingly S 
high for the structure, indicating, 
as Bergmann and Tschudnowsky 
conclude, that the rings do not lie S 
in the same plane. Indeed, a considerable folding 
of the structure about the two sulfurs must occur 
in order to give such a moment, for the C-S 
valences must act at such angles as partially to 
oppose one another’s moments. The reason for 
the discrepancies in the values for thianthrene 
and the mercaptans is presumably the difficulty 
of purification, but it is believed that the ma- 
terials finally used in the present investigation 
are reliable as to purity. Impurity should account 
for the difference found by Hunter and Partington 
between. the higher mercaptan and the ethy! 
compound, the value for which is nearer than 
the others to the values found in the present 
work for the higher compounds. 


°C. P. Smyth and W. S. Walls, J. Am. Chem. Soc. 53, 
527 (1931). 
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DIPOLE MOMENTS OF MERCAPTANS AND SULFIDES 


TaBLE I. Dielectric constants, densities and polarizations. 


€ d 
23° 25° 50° 25° 50° 


Benzene-thianthrene 


0.00000 2.276 2.226 0.8730 0.8461 (26.68 26.76=P,) 
.00624 2.301 2.249 .8782 .8514 106.8 103.7 
01051 2.318 2.264 .8817 .8549 106.6 103 8 
01053 2.318 2.264 .8816 .8549 107.4 103.7 
01199 2.324 2.270 .8828 .8561 106.8 104.3 
01351 2.330 2.275 .8840 .8572 107.4 104.5 
01701 2.342 2.288 .8868 .8600 105.5 104.4 
01761 2.344 2.289 .8872 .8605 105.6 103.4 
Benzene-ethy]l sulfide 
0.01815 2.334 2.278 0.8717 0.8449 79.6 76.9 
03496 2.389 2.326 .8711 .8443 79.0 76.3 
05603 2.458 2.386 .8700 .8432 78.8 75.7 
10292 2.615 2.523 .8677 .8410 78.2 75.2 
33000 3.358 3.183 .8572 .8310 73.5 71.7 
50070 3.946 3.705 .8498 .8239 71.5 70.0 
60390 4.313 4.002 .8457 .8197 70.2 68.7 
1.00000 5.723 5.236 .8312 .8056 66.3 65.5 
Benzene-n-amy] sulfide 
0.01269 2.315 2.261 0.8717 0.8451 107.9 105.6 
02397 2.349 2.291 .8705 .8441 108.0 105.2 
03627 2.384 2.322 .8692 .8430 107.7 104.8 
09009 2.531 2.453 .8643 .8395 106.8 103.9 
14166 2.665 2.568  .8606 .8362 106.7 103.5 
.28299 2.956 2.832 .8526 .8296 104.6 102.2 
46309 3.244 3.088 .8460 .8244 103.0 100.8 
1.00000 3.826 3.594 .8356 .8162 101.1 99.0 
Benzene-n-butyl mercaptan 
0.02771 2.353 2.294 0.8714 0.8446 72.9 70.4 
05204 2.421 2.355 .8703 .8436 72.4 70.2 
07853 2.498 2.423 .8686 .8421 72.7 70.3 
13773 2.666 2.571 .8659 .8396 71.6 69.4 
22957 2.925 2.803 .8617 .8358 70.2 68.3 
37172 3.315 3.142 .8562 .8304 67.9 66.2 
56844 3.847 3.611 .8492 .8239 65.4 64.1 
1.00000 4.952 4.586 .8368 .8123 61.2 60.4 


Benzene-n-amyl mercaptan 


0.02123 2.337 2.280 0.8718 0.8451 79.0 76.7 
02965 2.362 2.302 .8712 .8444 79.6 77.0 
05362 2.427 2.361 .8698 .8432 78.2 76.0 
06427 2.460 2.387 .8692 .8427 78.8 75.9 
13816 2.663 2.567  .8652 .8392 77.2 74.8 
21759 2.877 2.753 .8612 .8357 76.0 73.6 
23514 2.914 2.789 .8603 .8350 75.3 73.1 
40585 3.337 3.158  .8535 .8291 73.0 71.0 
51977 3.600 3.388  .8493 .8255 71.7 69.9 
1.00000 4.547 4.230 .8373 67.4 66.3 


TaBLE II. Refractions, orientation polarizations and dipole 


moments. 
P.—MRpb 
MRp 25° 50° 25° 50° 
Thianthrene 65.9 41 38 141 1.41 


Ethyl sulfide 
n-Amyl sulfide 
n-Butyl mercaptan 


n-Amyl mercaptan 


28.57 51.7 48.2 1.58 1.58 
56.44 51.9 48.8 1.58 1.59 
28.40 45.4 42.8 1.48 1.49 
33.03 47 44 1.50 1.51 
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TABLE III. Dipole moments of sulfides and mercaptans 
(108). 


Hydrogen sulfide 0.93* 


Mercaptan Sulfide 
W.andS. H.andP. W.andS. H. and P. 
Methyl 1.40 
Ethyl 1.39 1.58 1.57 
Propyl 1.33 135 
Butyl 1.48 1.32 1.57 
Amy! 1.50 1.58 


* J. B. Miles and C. T. Zahn, Phys. Rev. 32, 497 (1928). 


Hunter and Partington have based a fairly 
satisfactory calculation of the moments of the 
sulfides and the mercaptans upon the assumption 
of a 90° sulfur valence angle. A fairly satisfactory 
result is also obtained by using the 146° valence 
angle and the high mc-u1+mc-s moment, 2.5 
10-8, calculated by Smyth and Walls.'° The 
moment, 1.46X10-", thus calculated for the 
sulfides is very close to the experimentally 
determined values in Table III. However, as 
pointed out by Smyth and Walls, this large sulfur 
valence angle and surprisingly high C-S moment 
are dependent upon a somewhat inaccurate value 
in the literature for the moment of p-CIC,H,SC,H; 
and upon the assumption that the moment of the 
para-substituted group is approximately inde- 
pendent of the other moments in the molecule. 
The differences between the moments of para- 
substituted phenols, anisoles, etc.,'° make it 
appear probable that this assumption is not 
generally justifiable in the case of the ethers and 
thioethers, thus leaving the question of the sulfur 
valence angle open as far as dipole moments are 
concerned. 

Hunter and Partington have explained the 
differences between the moments of corresponding 
oxygen and sulfur compounds largely in terms of 
the considerable difference in polarizability of 
oxygen and sulfur, as was done roughly for water 
and hydrogen sulfide by Smyth. A C-S moment 
inherently larger than the H-S and a C-O 
inherently smaller than the H—-O, both of which 
probable facts are influenced by the sulfur and 
oxygen polarizabilities, can account for the 
differences-in moment. Not only the increase in 
moment from the hydrogen halides to the 


10C, P. Smyth and W. S. Walls, J. Am. Chem. Soc. 54, 
3230 (1932). 
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corresponding methyl halides, but, to a smaller 
extent, also the increase of the methyl halide to 
the corresponding ethyl compound increases with 
increasing polarizability of the halogen. It thus 
appears reasonable that, although the C—O 
moment is smaller than the H-—-O, the C-S 
moment should be larger than the H-S, which 
gives the mercaptans higher moments than H2S 
and the sulfides higher moments than the 
mercaptans. 

A surprising effect is apparent in the de- 
pendence of the observed polarization of the 
sulfide upon concentration. The highly as- 
sociated alcohols show a strong and complicated 
dependence of apparent polarization upon con- 
centration because of the strong effects of the 
H-O dipoles upon one another,'! while the ethers 
show a polarization independent of concentration, 
presumably because the alkyl groups screen the 
C-O dipoles in adjacent molecules from one 
another. As would be expected, the supposedly 
unassociated mercaptans show (see Table I) a 


1C, P. Smyth, Dielectric Constant and Molecular Struc- 
ture, p. 180, The Chemical Catalog Company, Inc., New 
York, 1931. 
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change of only 10 to 15 units of observed 
polarization between the pure liquid and the 
substance at infinite dilution, the difference from 
the behavior of the alcohols presumably being 
due to the partial screening of the H-S dipole by 
the structure of the sulfur atom. The sulfides in 
Table I, however, show a dependence of polari- 
zation upon concentration which, when the 
difference in molecular size is taken into account, 
is as great as that of the mercaptans. Thus, the 
ethyl sulfide molecule, which is approximately 
the same size as that of butyl mercaptan, in- 
creases in apparent polarization from the pure 
liquid to infinite dilution, 14.0 at 25° and 11.3 at 
50°, while butyl mercaptan increases 12.6 at 25° 
and 10.8 at 50°. It appears possible that the 
larger size of the sulfur atom as compared to the 
oxygen prevents complete screening of the C-S 
dipoles by the alkyl groups so as to permit a 
larger external field than in the case of the 
oxygen ethers and consequently some molecular 
interaction. The larger moments of the sulfides 
are also a factor, possibly, helping to make the 
intermolecular action of the sulfides as large as 
that of the mercaptans. 
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